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STABILITY RELATIONS OF SIDERITE (FeC0 3 ), DETERMINED 
IN CONTROLLED-f 0o ATMOSPHERES 


Bevan M French 
Planetology Branch 


ABSTRACT 

Stability relations of sidente(FeC0 3 )m the system Fe-C-O were determined between 500 and 
2000 bars m a CO a + CO atmosphere as a function of T, P F (= P C02 + P co ) and f 02 , using solid- 
phase oxygen buffers Siderite was synthesized for the experiments by decomposition of ferrous 
oxalate dihydrate (FeC 2 0 4 *2H 2 0)in sealed tubes at 2000 bars P H „ 0 and approximately 380°C Dur- 
ing stability runs, f 02 was controlled by surrounding the siderite sample with either graphite or 
hematite-magnetite buffer mixtures. 

The stability field of siderite + gas m P F - T - f 02 space is a narrow wedge bounded by stable 
divariant reaction surfaces representing decomposition of siderite to (1) hematite + gas, (2) 
magnetite + gas, (3) magnetite + graphite + 0 2 . The latter reaction surface lies entirely in the 
condensed region below the graphite buffer and is not experimentally accessible. Reactions of sid- 
erite to form iron or wustite (Fe 1 _ x O) are not stable 

Points along the umvanant equilibrium curve siderite + hematite + magnetite + gas (SHMG) 
were located reversibly at P F = 500 bars, T = 363 ± 10°C, log f o 2 = -24.8, P F = 1000 bars, 
T = 365 ± 10°C,logf 0 = -24.7 The umvanant equilibrium curve siderite + magnetite + graphite 
+ gas (SMGrG)passes through the following points P F = 500 bars, T = 455 ± 10°C, logf 0 = -25.7; 
P F = 1000 bars, T = 458 ± 10°C, logf o 2 = -25 2, P F = 2000 bars,T = 465 ± 10°C,logf 0 , = -24.4. 

hi an isobaric f 0 ,-T section, the umvanant curve siderite + hematite + gas (SHG) corresponds 
to low temperatures and relatively high f 0 values The umvanant curve siderite + magnetite 
+ gas (SMG) is stable at higher temperatures and relatively lower f 02 values. The stability of the 
assemblage siderite + gas is strongly dependent on T and f 0a , but is practically independent of the 
value of P C q 2 + P co* The umvanant equilibrium siderite + magnetite + graphite + gas fixes the 
maximum stability temperatures of siderite, 455-465°C in the range P F = 500-2000 bars. In this 
range, siderite is not stable at f 02 values above about 10~ 24 bars and decomposes to either hematite 
(lower T) or magnetite (higher T). The determined stability field of siderite + gas lies at f 02 values 
above the stability field of fayalite (Fe 2 Si0 4 ) and therefore reactions between siderite and quartz 
to form fayalite are not stable 

The stability relations of natural iron-rich carbonates are more complex than those deter- 
mined for pure siderite. The presence of other volatile species in a natural gas phase, producmg 
a condition in which P co , + P c0 < P F , will produce slightly lower decomposition temperatures, and 
the presence of H 2 Q may stabilize iron hydroxide phases at lower temperatures. By contrast, the 
substitution of Mg 2+ , Mn 2+ , or Ca 2+ for Fe 2+ will stabilize iron- rich carbonates at higher tem- 
peratures and f 02 values. 

During metamorphism, natural iron-rich carbonates react with available quartz and water to 
produce iron-rich amphiboles (grunerite). The present experimental data, combined with geological 
studies of metamorphosed iron-formations, suggest that this reaction occurs at about 300-450°C. 
Conversely, the coexistence of siderite + quartz m many hydrothermal veins indicates tempera- 
tures below this value. Conversion of quartz-free siderite units to magnetite in contact- 
metamorphie aureoles requires temperatures exceeding 400-450°C. 
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SYMBOLS AND NOTATIONS USED 

Symbols used generally follow conventions of Thompson (1955) and Greenwood (1961). 

T = temperature (°K). 

P s = isotropic pressure on solid phases 
P F = total pressure of the fluid (gas) phase. 

P E = "osmotic equilibrium" pressure of component 1 m a fluid phase (Greenwood, 1961). 

. * 

P t = total pressure in a closed experimental system where P t = P F = P s . 

P, = partial pressure of component i in a closed experimental system where 2P X = P t . 
f 1 = fugacity of component i m a fluid phase. 

y 1 = fugacity coefficient of component i. in the fluid phase, defined by ■y i = f 1 /P 1 . 

AV s = volume change of the solid phases in a reaction, defined as AV s = 2(v s ) product; . 

1 reactants 

AG f ° T = standard Gibbs free energy of formation of a compound from its elements at a speci- 
fied temperature and unit fugacities of gaseous components 

AG X ° = standard Gibbs free energy of reaction, equal to 2(AG f ° x ) products - 2(AG f ° T ) reactaats . 

AH f ° T = standard enthalpy of formation of a compound from its elements. 

= standard enthalpy change in a reaction. 

K x (T) = equilibrium constant of reaction i , written with fugacities of gaseous components, a 
function of T only. 

c p = molar heat capacity of a substance at constant pressure 
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ABBREVIATIONS FOR PHASES IN TABLES AND ILLUSTRATIONS 


S 

H 

'M 

I 

G 

G* 


= sidente (FeC0 3 ) 

= hematite (Fe 2 0 3 ) 

= magnetite (Fe 3 0 4 ) 

= iron (Fe) 

= gas phase m a C0 2 + CO atmosphere 

= gas phase (0 2 ) in condensed region below the graphite buffer curve. 


Wu = wustite (Fe^O) 
Gr = graphite (C) 

Q = quartz (Si0 2 ) 

F = fayalite (Fe 2 Si0 4 ) 


hi tables of experimental results, occurrence of trace amounts of a phase are indicated by 
parentheses, as (S). Such occurrences are not believed to represent stable reaction products and 
may result from contamination or from metastable reaction. Clearly metastable products m a 
significant amount are indicated by an asterisk, as (S*). Small letters (m, s, h, etc.) designate 
phases which may be in equilibrium with the observed phases but which were not detected, e.g , the 
assemblage S + M + h. 
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STABILITY RELATIONS OF SIDERITE (FeC0 3 ), DETERMINED 

IN CONTROLLED-! 0 ATMOSPHERES 

2 


INTRODUCTION 

Iron-bearing minerals are of particular interest in petrogenetic studies, because they may 
participate in reactions involving oxidation and reduction (Eugster, 1959, Buddmgton and Lmdsley, 
1964). The occurrence and composition of iron-bearing minerals thus reflect in part the value of 
f o in the environment during their formation. Knowledge about the stability of iron-bearing min- 
erals is thus directly applicable to analyses of conditions of formation of iron-bearing assemblages 
and to discussions of the behavior of volatile components (H 2 0, CO,, 0 2 et al.) during metamorphism 
(Mueller, 1960, Kranck, 1961, French, 1968, Klein, 1966, Butler, 1969). 

Siderite (FeCO a ) is a common constituent of low-grade sedimentary iron formations and of 
some hydrothermal veins, and information about its stability Is of great importance in discovering 
how such deposits have formed and what were the compositions of the volatile phases with which 
they coexisted Unlike most other carbonates, the stability of siderite is dependent on temperature, 
pressure, and the fugacities of both C0 2 and 0 2 , and exact determination of its stability requires 
simultaneous control of both f co and f 0 . 

Two types of decarbonatxon reactions exist for natural carbonate minerals (Jamieson and 
Goldsmith, 1959) 


RC0 3 = RO + C0 2 , (1) 

\ 

aRC0 3 +|o 2 = R a O a+1 + aC0 2 (2) 

In reaction (1), the divalent cation (R) is not oxidized, and equilibrium is independent of the partial 
pressure of oxygen In reaction (2), the cation may assume different valences and form a variety 
of oxides at different values ol Po 2 . The stability of a carbonate which decomposes by reactions 
of the second type will depend on the partial pressure of oxygen as well as on temperature and P C o 2 • 

Most experimental studies have dealt with reactions of the first type Stability relations as a 
function of P cc , 2 and T have been determined for calcite (CaC0 3 ) and magnesite (MgC0 3 ) (Harker 
and Tuttle, 1955a) and for smithsomte (ZnC0 3 ) (Harker and Hutta, 1956). The analogous reaction 
of dolomite to form calcite, periclase (MgO), and CO a has been determined by Harker and Tuttle 
(1955a) 

Natural carbonates which decompose according to reaction (2) are siderite (FeC0 3 ) and rho- 
dochrosite (MnC0 3 ). Their stability is thus a function of T, P COa , and Pq 2 , and these three vari- 
ables must be simultaneously controlled for rigorous stability determinations Such control has 
been made possible by the use of solid-phase oxygen buffers m an atmosphere of C0 2 + CO (French, 
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1964a, French and Eugster, 1965) Preliminary data have been presented for the stability of sid- 
erite (French, 1965, French and Eugster, 1965, French and Rosenberg, 1965) and a detailed study 
of rhodochrosite stability has recently appeared (Huebner, 1969). Siderite decomposes stably to 
hematite and magnetite by oxidation reactions of the second type. The rhodochrosite stability field 
is bounded chiefly by analogous oxidation reactions, but rhodochrosite also decomposes stably to 
manganosite (MnO) over a limited range of (Huebner, 1969) 


DESCRIPTION OF THE CRYSTALLINE PHASES 

The term siderite refers here, to calcite-type carbonates in which the FeC0 3 end-member is 
dominant. In the carbonate crystal structure, complete solid solution apparently exists between 
Fe 2+ and Mg 2+ and between Fe a+ and Mn 2+ (Palache et al , 1944, Goldsmith, 1959, Rosenberg, 1963a, 
1963b, 1967) A wide miscibility gap exists, however, between CaC0 3 and FeC0 3 . Limited solu- 
bility of Ca 2+ an siderite is observed, the- experimentally determmed solubility does not exceed 10 
mole percent CaC0 3 at temperatures as high as 550°C (Rosenberg and Barker, 1956, Rosenberg, 
1963a, 1963b, 1967). Fe 2 * is more soluble in calcite, at 550°C, caleite may contain 20 -mole percent 
FeC0 3 (Rosenberg, 1963a). 

Although carbonates with over 95 mole percent FeC0 3 have been reported (Hutchinson, 1903, 
Ford, 1917), most natural siderites contain significant amounts of Mn 2 ' and Mg 2 *. Analyses of 
natural siderites and related carbonates have been compiled by Ford (1917), Sundius (1925a, 1925b), 
Schokhtsch (1935), Palache et al. (1944), Hugi (1945), Smythe and Dunham (1947) and Deer et al. 

(1962, p. 272-277). The composition of siderite m the Lake Superior iron formations varies be- 
tween 74 and 92 mole percent FeC0 3 (James, 1954), and analyzed siderites from unmetamorphosed 
iron formation on the Mesabi Range contain between 61 and 80 mole percent FeCO a , 18-31 mole 
percent MgC0 3 , and 2-7 mole percent CaC0 3 (French, 1968). 

V. 

In natural occurrences, siderite is commonly associated with iron-rich dolomite (ankente). 
Siderite is most commonly found m sedimentary iron formations (James, 1954, Pettijohn, 1957, 
Carrozi, 1960) and in hydrothermal veins (Lindgren, 1933, Wahlstrom, 1935, Legraye, 1938, Fabian 
et al , 1957; Shaw, 1959) The dominant carbonates in carbonatite bodies associated with alkalic 
igneous rocks are generally iron-poor, but both siderite and ankente have been reported from 
several localities (Pecora, 1956, Smith, 1956). Siderite has also been observed in salt-dome cap 
rock (Rolshausen, 1934), in a bauxite-bearing clay (Burchard, 1924), and in a tourmaline -bearing 
pegmatite (Simpson, 1929). Siderite is present in moderately metamorphosed iron-rich rocks, 
generally associated with the iron-amphibole grunente (Tilley, 1938, Gustafson, 1933, James, 1955, 
Klein, 1966), but has not been reported from highly-metamorphosed rocks containing iron-rich 
pyroxenes (James, 1955, Gundersen and Schwartz, 1962, French, 1968, Bonnichsen, 1968). 

- Magnetite, Fe 3 0 4 , is the most commonly observed member of the spinel group. Several ca- 
tions, particularly Ti, Cr, Al, and Mg, may substitute for Fe m the structure (Palache et al., 1944.). 
Above 900°C, magnetite may depart from stoichiometry by cation deficiency (Greig et al., 1935, 
Darken and Gurry, 1945, 1946). Magnetite is a common accessory mineral m sedimentary, igneous. 
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and metamorphic rocks, in hydrothermal veins, and m contact-metasomatic deposits Magnetite 
has also been recognized as an original constituent of sedimentary iron formations (James, 1954), 
and is generally present when such rocks are metamorphosed (LaBerge, 1964, French, 1968) 

Hematite, a-Fe 2 0 3 , is a common mineral in igneous and metamorphic rocks * It is the dominant 
iron oxide m post-Precambrian sedimentary iron formations Hematite can contain a small amount 
of A1 2 0 3 m solid solution (Turnock and Eugster, 1962); Ti, Mn, and other cations may also be pres- 
ent Other polymorphs of Fe 2 0 3 (see Rooksby, 1961) were not observed in the present study 

Graphite, the hexagonal crystalline form of carbon, is found in igneous rocks, m metamorphosed 
sediments and coals, and m vein deposits of uncertain origin Much carbonaceous material from 
unmetamorphosed sediments, which has been considered graphite, is probably amorphous or poorly- 
crystalline carbon or hydrocarbons A thorough summary of the properties of graphite and related, 
compounds is given by Ubbelohde and Lewis (1960). Graphite may form during metamorphism from 
noncrystallme or poorly crystalline carbonaceous matter in sediments (Quinn and Glass, 1958, 
French, 1964b) 

PREVIOUS STUDIES OF SIDERITE STABILITY 

Instability Under Atmospheric Conditions 

The instability of siderite under the conditions of low P cc>2 and high P D present at the sur- 
face has been demonstrated by (1) thermodynamic calculations, (2) observed alteration of natural 
siderite, (3) experimental oxidation of siderite at room temperature 

Thermodynamic calculations of the stability of siderite, hematite, and magnetite indicate that 
siderite is not stable under the values of P 0 (0.2 atm) characteristic of the atmosphere (Holland, 
1959^4.965, Garrels, I960) At the value of P co in the atmosphere (3 x 10 -4 atm), and at 25°C, 
siderite should decompose to hematite above 10" 75 atm P^ 

Natural siderite, exposed to the atmosphere or to ground water, commonly oxidizes to hematite, 
limonite, or goetlkte, and more rarely to magnetite (Palache et al , 1944, Petti] ohn, 1957, p 145, 
Burchard, 1924, Smythe and Dunham, 1947, Deer et al., 1962)' Silliman (1820) described the oxida- 
tion of pieces of siderite exposed for about 50 years on the mine dump at Roxbury, Connecticut, a 
rim of dark, oxidized material one-quarter tc^ one-half inch thick surrounded a core of lighter un- 
altered siderite 

Siderite oxidizes readily m laboratory experiments at room temperature A siderite sample 
stored at room temperature underwent progressive oxidation (Schaller and Vlisidis, 1959) over a 
period of 43 years, the FeO content decreased linearly from 59 42 to 0 74 percent. Jamieson and 
Goldsmith (1960) oxidized natural siderite to hematite by prolonged grinding at room temperature 
Synthetic siderite partially oxidized in air during drying at 100°C (Sharp, 1960) A kmetic study of 
siderite oxidation (Segum, 1966) between 25°C and 300°C suggests that the process is a first-order 
reaction with an activation energy of about +10 kcal/mole 

DTA Studies 

Previous experimental studies of the decomposition of siderite have been carried out 
largely by differential thermal analysis (DTA) or by thermo- gravimetric analysis (TGA). The 


3 



DTA pattern of sidente obtained by many investigators (Cuthbert and Rowland, 1947, Kerr and 
Kulp, 1947, Fredenekson, 1948, Rowland and Jonas, 1949, Beck, 1950, Kulp et al , 1951, Kissinger 
et al., 1956, Powell, 1965) consists of an exothermic peak in the range 400-650°C, followed by an 
endothermic peak at 750~850°C The first peak has been ascribed to the breakdown of sidente to 
FeO, the second to the oxidation of FeO to magnetite or hematite This interpretation is uncertain 
because stoichiometric FeO does not exist and the compound Fe,_ x O (wustite) is stable only above 
560°C and at low values of P 0 ^ {Darken and Gurry, 1945, 1946) 

The DTA studies have been made in atmospheres in which P 0 is unknown and is probably too 
high for equilibrium decomposition of siderite The exothermic peak may indeed represent the 
non- equilibrium reaction of sidente to an iron oxide (hematite, magnetite, or possibly wustite), 
while the higher-temperature endothermic peak may indicate complete conversion of the iron oxides 
to hematite (Powell, 1965) 

The interpretation of the DTA results is uncertain. The partial pressure of oxygen was not 
controlled and is not known, it was probably much higher than the values required for equilibrium 
between siderite and the iron oxides. This interpretation is favored by the observed strong effects 
of different furnace atmospheres on reaction temperatures (Rowland and Jonas, 1949, Kissinger 
et al., 1956, Powell, 1965) and by the observed lowering of the temperature of the exothermic peak 
with decreasing heating rate (Kissinger et al , 1956) 

Hydrothermal Studies 

The stability relations of side rite-bearing carbonate assemblages in the subsolidus region 
of the system CaC0 3 -MgC0 3 -FeC0 3 -MnC0 3 have been studied experimentally at temperatures 
up to 550°C and fluid pressures as high as 4 kb (Rosenberg and Harker, 1956, Rosenberg, 
1963a, 1963b, 1967). Similar assemblages in the subsolidus region of the system CaC0 3 - 
MgC0 3 -FeC0 3 have been examined at temperatures up to 800°C at total pressures of 15 kb (Gold- 
smith et al., 1962). In these studies, the value of P C c 2 was maintained high enough to prevent dis- 
sociation of the carbonates. The value of was not strictly controlled and significant alteration 
of siderite to magnetite was occasionally observed (Rosenberg, 1960, p 8, Goldsmith et al., 1962, 

p. 660). 

More recent experimental studies have concentrated on the stability and decomposition of 
siderite itself Segum (1968a), in a study of the complicated Fe-C-0-S-H 2 0 system, decomposed 
siderite to iron oxides by various techniques, both with and without the use of solid-phase oxygen 
buffers. In buffered experiments, he obtained decomposition temperatures at a total pressure of 
2 kb which ranged from 296°C (HM buffer) to 438°C (MI buffer) His higher temperatures are clearly 
metastable with respect to the graphite buffer (French and Eugster, 1965), and interpretation 
of his other results is difficult because of the complexity of the system in which he worked and the 
lack of demonstrated reversibility m his reactions. 

Johannes (1968, 1969) investigated the stability of Fe-Mg carbonates in equilibrium with Fe 2 + - 
and Mg 2f -bearing aqueous solutions While the value of was not explicity controlled m his 
experiments, he observed that siderite and Fe-Mg carbonates containing more than 30 mole per- 
cent FeCQ 3 decomposed to magnetite at about JOO^C 
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A study of siderite decomposition relations in the system Fe-C-O (Weidner and Tuttle, 1964, 
Weidner, 1968) was carried out using sealed tubes m a CO a atmosphere at pressures up to 10 kb. 
Temperatures determined for the equilibria siderite + hematite + magnetite + gas and siderite 
+ magnetite + graphite + gas are higher than those determined m the present study (French, 1965, 
French and Rosenberg, 1965), and a break m the slope of the curve for the equilibrium siderite 
+ magnetite + graphite + gas, located at 5 kb, was interpreted as indicating a polymorphic transi- 
tion of siderite to a high-pressure phase (Weidner and Tuttle, 1964). 

The present investigation of siderite stability is the first in which explicit control of the value 
of P 0j is obtained by use of open capsules combined with solid-phase oxygen buffers (French and 
Eugster, 1965, Huebner, 1969). A detailed discussion of the difference in results obtained on sid- 
erite stability by different experimental methods of different investigators is given below (p 31). 


THERMODYNAMIC CALCULATION OF SIDERITE STABILITY 

Sufficient thermodynamic data exist for calculation of the stability field of siderite m terms of 
P co , P 0 , and T (see, e g , Garrels and Christ, 1965, Holland, 1959, 1965, Mel'nik, 1964, Yui, 1966). 
The stability field of siderite is considered to be bounded by the following possible reactions 


2FeC0 3 + 

2 °! ~ Fe 2°3 + 2C0 2 


siderite 

= hematite + gas 


log K 3 (T) 

= 2 lo g f co 2 " \ lo e f o 2 

(3) 

3 FeC0 3 + 

I°2 = Fe 3 0 4 + 3CO, 


siderite 

= magnetite + gas 


log K 4 (T) 

= 3 log f C0 2 - \ log f 0 2 

(4) 


0 947 FeC0 3 + 0 026 0 2 = Fe Q 9+? 0 + 0 947 C0 2 

siderite - wustite + gas 

log K s (T) = 0 947 logf^ - 0 026 logf^ - (5) 

1 

FeC0 3 = Fe + C0 2 + 0 2 

siderite = iron + gas 

logK 6 (T) = log f C o 2 + \ log f ^ (6) 
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( 7 ) 


5 

3 FeCOg = Fe 3 0 4 + 3C + tj- 0 2 

siderite = magnetite + graphite + gas 
logK 7 (T) = |logfo 2 
3 

FeCOg = Fe + C + 0 2 

siderite = iron + graphite + gas 
logK 8 (T) = -| lo g f o 2 

Reactions (7) and (8) can be stable only at f values below that of the graphite buffer, they lie 
entirely within the condensed region and do not involve C0 2 or CO (French and Eugster, 1965). Cal- 
culated stability fields of siderite, bounded by reactions (3), (4), and (6), have been presented by 
Garrels (1960) for 298 °K and by Holland (1959) for 300-600°K. Neither author considered possible 
reactions involving siderite and graphite. 

In the present study, AG f ° T data for siderite were calculated from thermodynamic data (Kelley 
and Anderson, 1935). Values of AG t ° for each reaction were calculated using values for the oxides 
tabulated by Coughlin (1954) The calculated values of ag t ° are essentially linear functions of T from 
298° to 1200°K (Table 1). Uncertainties in the values were calculated from the stated uncertainties 
of Coughlin (1954) for the oxide data, and a 1 percent uncertainty was assumed for the siderite data, 
the total uncertainty in AG X ° is thus only a minimum value. 

Table 1 


Calculated Thermodynamic Data for Siderite Decomposition Reactions 


Reaction 

AG 2 ° 9S (cal ) 

AG t c {cal ) 

logK, <T) 

(3) siderite + O a = hematite + C0 2 
2 FeC0 3 + 1/2 0 2 = Fe 2 0 3 + 2 C0 2 

'-43,850 ±4,500 

-26,500 - 54 29 T 

+6012/T + 11 87 

(4) siderite + 0 2 = magnetite + C0 2 
3 FeC0 3 + 1/2 0 2 = Fe 3 0 4 + 3 C0 2 

-42,870 =b 6,200 

-13,900 - 96 44 T 

+3039/T + 21 08 

(5) siderite + 0 2 = wustite + C0 2 

0 947 FeCOg + 0 026 0 2 = Fe 0 947 0 + 3 C0 2 

+4,490 ± 1,800 

+16,300 - 39 6 T 

-3564/T + 8 66 

(6) siderite = iron + C0 2 + 0 2 
FeCOg = Fe + C0 2 + 1/2 0 2 

+66,780 ± 1,700 

+84,600 - 59 65 T 

-18,495/T + 13 04 

(7) siderite = magnetite + graphite + 0 2 
3 FeCOg = Fe 3 0 4 + 3 C + 5/2 0 2 

+239,910 ± 6,100 

+268,300 - 91 69 T 

-58,656/T + 20 05 

(8) siderite = iron + graphite + 0 2 
FeCOg = Fe + C + 3/2 0 2 

+161,040 ± 1,600 

+178,700 - 58 97 T 

-39,0G8/T + 12 89 
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Values of as° and AH 0 for each reaction (Table 1) were calculated from the free energy data, 
using the relations ~ 


/ <9AG°\ 

( <JT J p = “*5° 

(9) 

/ AG°\ 


1 q— | = AH° 

(10) 

\ "t / 

\ 1 /p 



(Lewis et al., 1962, p 165), AS° and AH° for all reactions axe virtually independent of T over the 
temperature range 2S8-1200°K 

The equilibrium constants K(T) for each reaction were calculated from the relation 

AG t ° = -2 303 KT log K(T) 


logK(T) 


~AG° 

2 303 RT 


(ID 


and are given m Table 1 

Values of equilibrium f 0 as a function of T can then be calculated from the equilibrium constant 
for any specified total pressure. In those reactions in which C0 2 also participates, it is reasonable 
to assume that P F = P c0 for temperatures below 600°C and f 0 values above those of the graphite 
buffer (French and Eugster, 1965). Values of f co were calculated for these reactions from the 

2 i 

relation f 1 = y i P t , using fugacity coefficient data tabulated by Robie (1962). The results of such 
calculations for P F = 2000 b are shown in Figure 1 with the curves for several specified oxygen 
buffer assemblages 

The calculated stability field of siderite thus obtained is a narrow wedge bounded by reac- 
tions (3), (4), and (7), with stable decompositions to hematite, to magnetite, and (below the graphite 
buffer curve) to magnetite + graphite. In a C0 2 + CO atmosphere, the breakdown of siderite to iron 
(reaction 6) lies m the condensed area below the graphite buffer curve where C0 2 and CO do not 
exist (French and Eugster, 1965). This reaction is therefore metastable, and the siderite stability 
field is thus much smaller than those presented by Holland (1959) and Garrels (1960). The break- 
down of siderite to wushte is also metastable, m agreement with the calculations of Muan (1958) 
for 1 atm total pressure The stability N field of siderite does not intersect that of fayalite (Fe 2 Si0 4 ), 
the upper limit of which is bounded by the quartz + fayalite + magnetite (QFM) buffer curve (Fig- 
ure 1). Reaction between siderite and quartz to produce fayalite is not stable at these pressures 
(French and Eugster, 1962, French, 1964a). 

The system Fe-C-O contains three components, from the Phase Rule, F = 5 - P Reactions (3) 
and (4), in which two solid phases are m equilibrium with a gas phase, are therefore divariant 
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Temperature, °C 

Figure 1— Isobanc section at P F = 2000 bars, showing sidertte stability relations cal- 
culated from thermodynamic data The stable field of siderite (narrow ruling) is 
bounded by decomposition reactions to hematite (SHG), magnetite (SMG), and mag- 
netite + graphite (SMGrG - *). The latter reaction lies entrrely in the condensed re- 
gion below the graphite buffer (Gr) Because of the reactions involving graphite, this 
field is much smaller than that calculated using reactions not involving graphite 
(Holland, 1959, Garrels, 1960) Other solid lines (HM, Gr, QFM, and Ml) designate 
solid-phase oxygen buffers 

reaction surfaces in P„ - f n - T space Reaction (7), in which three solid phases and a gas co- 

2 

exist, is umvariant and lies entirely within the region of the diagram where P F = P 0 (French and 
Eugster, 1965), 

We may imagine addition of an inert ideal gas to the condensed region in such a manner that 
the value of P F below the graphite buffer curve is always kept equal to the value of P F = P co ^ + P co 
above it Under these conditions, the region below the graphite buffer curve is a four- component 
system Fe-C-O-gas and the four-phase assemblagebecomesadivariant surface m P F - P 0 - T 
space. The surface will have an orientation similar to that of other buffering assemblages (Eugster 
and Wones, 1962). 

With this reservation, reaction (7) likewise becomes a divaraant reaction surface m P F - P 0 -T 
space Tiie three divariant surfaces of reactions (3), (4), and (7) intersect m two umvariant curves 
The curve siderite + hematite + magnetite + gas (SHMG) is fixed by the intersection of reac- 
tions (3) and (4) with the hematite- magnetite (HM) buffer The curve siderite + magnetite + 
graphite + gas (SMGrG) is the mutual mtersection of reactions (4) and (7) with the graphite buffer _ 
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surface This latter univariant curve is the only part of the surface of reaction (7) which is experi- 
mentally accessible These univariant curves appear m Figure 1 as isobaric invariant points, 
SHMG at 228°C and SMGrG at 283 °C. . (A third intersection, corresponding to a univariant curve for 
the equilibrium siderite + hematite + graphite + gas, is theoretically possible, but appears to be 
metastable under reasonable experimental conditions.) 

In the system Fe-C-O, five phases must coexist at a true invariant point (four solids plus gas) 
The point siderite + magnetite + wustite + iron + gas is apparently metastable, since it requires a 
minimum temperature of 560°C It is uncertain if the equilibrium siderite + hematite + magnetite 
+ graphite + gas is metastable or whether it is stable at temperatures too low to be meaningful or 
experimentally accessible. 

These thermodynamic calculations are useful for delineating stability relations m advance of 
actual experimental work However, because of the uncertainties inherent in the basic thermody- 
namic data, these calculated stability diagrams can be only approximations The relative stability 
relations can generally serve to indicate closely the stability field determined experimentally, but 
great differences between calculated and experimentally determined equilibrium temperatures can 
be expected. The effect of such uncertainties is indicated m Figure 2, a plot of AG against T for 
reactions (3) and (4) at total pressures of 1 bar and 2 kb. Values of ag are calculated from the 
relation 


AG = AG° 1- 2 303 RT^a logf co ^ +\ log (12) 

(a is the coefficient of siderite m each reaction, a - 2 for reaction (3) and a = 3 for reaction (4)). 
Values of log f are those for the hematite -magnetite (HM) buffer. Ideally, the two curves for 
reactions (3) and (4) should intersect each other and the line AG = 0 at the equilibrium temperature 
of the isobaric invariant point siderite + hematite + magnetite + gas The estimated uncertainty m 
AG" is the sum of the uncertainties mAG f ° T for the individual phases, the values are ±4 5 kcal/mole 
for reaction (3) and ±6 0 kcal/mole for reaction (4). These uncertainties produce an uncertainty in 
the equilibrium temperature of ±50°C-at P F = 1 bar and ±90°C at P F = 2 kb. Even these large un- 
certainties in equilibrium temperature are only minimum values. 

EXPERIMENTAL METHOD 
Introduction 

In experimental studies at low pressure, control of f 0 in a gas phase can be attained by using 
a gas mixture with a specified H 2 0/H 2 or C0 2 /C0 ratio that corresponds to the desired value of 
f 0 . This gas-mixing method is difficult to use at the high pressures characteristic of many de- 
hydration and decarbonation reactions of geologic interest 

In high-pressure hydrothermal studies of oxidation-reduction reactions, control of f 0 is there- 
fore attained by surrounding the sample with an "oxygen buffer" consisting of one or more solid 
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Temperature, °C 



Figure 2™Plot ofAGvs T for the reactions siderite = hematite + gas (5HG) and 
siderite - magnetite + gas (SMG) at total pressures of 1 bar and 2000 bars, at values 
of log f 0 ^ fixed by the hematite-magnetite buffer The calculated equilibrium tem- 
perature is defined by the intersection of the two curves with each other and with 
the AG = 0 line. Note that the uncertainties in AG (stippled areas) produce an area 
of uncertainty at the intersection (diamond-shaped ruled area) which corresponds to 
a large uncertainty in the equilibrium temperature Uncertainties are ±4 5 kcal for 
the SHG reaction and ±6 2 kcal for the SMG reaction (Table 1) 
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phases whose presence fixes' the oxygen fugaeity for any given temperature and total pressure 
(Eugster, 1957, 1959, French and Eugster, 1965). This technique has been successfully applied to 
hydrothermal studies of iron-bearing minerals whose stability is dependent on f 0 iron-spinels 
and -chlorites (Turnock, I960, Turnock and Eugster, 1962), iron-amphiboles (Ernst, 1960, 1961, 
1962), and iron-biotites (Eugster and Wones, 1962, Wones and Eugster, 1965) hi these studies, 
equilibrium m the gas phase was attained by diffusion of hydrogen through sealed platinum capsules 
which contain a sample surrounded by the buffer The attainment of such equilibnum was demon- 
strated by the consistence and reversibility of the reactions studied (Eugster and Wones, 1962, 
Turnock and Eugster, 1962), 

The buffering principle is quite general and the techniques can be applied to any system in- 
evolving solid- gas equilibria. Similar buffers have already been developed for use m C0 2 atmos- 
pheres (French and Eugster, 1965) and buffers controlling CH 4 and F'have been used experimentally 
(Eugster and Skippen, 1968).. 

The present method for studying the stability of siderite is based on the use of a C0 2 + CO at- 
mosphere, combined with solid polyphase oxygen buffers to control the value of f o 2 (French and 
Eugster, 1965, Huebner, 1969). Independent control of the two volatiles C0 2 and 0 2 is thus possible, 
since, for a given P F = P co + P co , the values of f co , f co , and f 0 for any given temperature can 
be calculated (French, 1964a, French and Eugster, 1965). For the range of temperature and f 0 
values investigated m the present study, If co is virtually equal to P F (French, 1964a). 

In this study, values of f 0 specified by the hematite-magnetite (HM) buffer are those calculated 
by Eugster and Wones (1962, p. 90), according to the equation 

log! = - 24, 912/T 1 14 41 + 0 019(P - 1)/T (13) 

U 2 t 

A slightly different equation, obtained more recently for the hematite-magnetite' buffer for tempera- 
tures above 750°C (Haas, 1968) produces no significant changes in f 0 values calculated for the 
present experiments. Data for the graphite buffer have been calculated previously (French and 
Eugster, 1965). For values of P F > 500 bars, the following approximation is satisfactory 

iog f = - 20, S86/T - 0 044 + logP F - 0 028(P t - 1)/T (14) 

Experimental Procedure 

The powdered samples, generally siderite, were contained in a small (3 mm O.D.) silver tube 
(Figure 3), this inner tube was surrounded by a buffer mixture contained m a larger silver tube 

i 

(4 5 mm O.D ) The ends of both tubes were crimped but not sealed. During a run, C0 2 from the 
surrounding atmosphere-diffused through the oxygen buffer, the CO a /CO ratio of the gas was thereby 
adjusted to the equilibrium f 0 specified by the buffer and the run temperature. 

Runs were made in Tuttle-type cold-seal "test-tube" type bombs (Tuttle, 1948, 1949) made of 
Haynes Stellite 25 alloy The bombs were externally heated m horizontal furnaces. Liquid CO 2 
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to pump 



Figure 3— Cutaway sketch of run assembly In most 
runs, the sample (large dotted area) was synthetic sid- 
erite, and the buffer (horizontal ruled area) was a 
mixture of either hematife+ magnetite or of magnetite 
+ graphite. Buffering of the sample is established by 
diffusion of the gas phase (stippled area) though the 
buffer. 


under pressure was supplied to the bomb from a 
commercial siphon-equipped tank by a modified 
Sprague pump operated by compressed air A 
maximum C0 2 pressure of about 30,000 psi can be 
produced by the pump, operating on 80-100 psi of 
compressed air. During the early parts of this 
study, pressure was measured with several gauges 
calibrated against a 14-inch Heise gauge Later, 
these were replaced by a self- calibrating Foxboro 
gauge. Recorded pressures are believed accurate 
to ±5 percent. 


Temperatures were measured by chromel- 
alumel thermocouples j[Flgure 3), using a Leeds 
and Northrop potentiometer. Each furnace -bomb 
combination was calibrated with National Bureau 
of Standards samples of tin (M P 231.9°C) and 
zinc (M P. 419 5°C), consistmg of metal filings 
vacuum-sealed in 5,5 mm O.D. Pyxex tubes The 
calibration runs indicate that the temperature er- 
ror is less than ±3°C, variation in run tempera- 
ture during a two-week run generally did not ex- 
ceed ±4°C, The maximum temperature uncertainty 
# is thus ± f 7°C, and the run temperatures are believed 
accurate to ±5°C 

Analytical grade Fe 3 0 4 and Fe 2 0 3 were used 
for buffers. Graphite was prepared by grinding 
spectrograph electrodes. The siderite used as 
samples was synthesized as needed (see below, 
p 13) and used soon after synthesis to avoid pos- 
sible oxidation during long storage at room 
temperature. 

At the end of a run, the bomb was isolated 
from the line and quenched m water. Pressure 
and temperature fell together during the quench 
period of a few minutes. Since runs of about two 
weeks duration were required to produce detec- 
table reaction in the charge, it is believed that no 
significant reaction occurred during thfe heating 
and cooling periods 
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In the majority of runs, the samples consisted of synthetic siderite, while the buffers were 
either mixtures of hematite and magnetite (HM) or of magnetite and graphite (MGr). In the latter 
case, the magnetite is not part of the buffer assemblage, because buffering is produced by the 
graphite alone (French and Eugster, 1965), but this arrangement allowed reversibility of the reac- 
tion to be demonstrated in a single run. If the run temperature was above the stability field of 
siderite, the siderite sample decomposed to the oxides (Plate IB) while the buffer remained un- 
changed. Within the stability field of siderite, the sample remained unchanged while siderite 
formed from the oxides in the buffer. The addition of magnetite to the graphite buffer allowed 
this reverse reaction to be detected. 

Because the sample is surrounded during the run by potential reaction products in the buffer, 
care was taken to avoid contamination. After the run, the sample tube was carefully extracted and 
cleaned; the ends of the sample tube were clipped and discarded. A sample of the buffer adjacent 
to the sample tube was examined separately from the remainder of the buffer. 

Sample and buffer were examined under a binocular microscope, then studied by X-ray dif- 
fraction and oil immersion methods. In many of the runs, particularly those with the HM buffer, 
the amount of reaction was slight. Reaction of less than 5 percent could not be detected by X-ray 
diffraction. Sample decomposition in small amounts could be detected under the binocular micro- 
scope by the red or gray-black color produced by the iron oxides and by the magnetism of the 
sample. Under the petrographic microscope, small amounts of sample breakdown was indicated by 
the presence of strongly magnetic grains or of larger grains of siderite showing discoloration or 
opaque inclusions (Plate IB). Buffer decomposition (growth of siderite) was indicated by the forma- 
tion of white areas in the buffer and could be verified under the petrographic microscope by the 
presence of birefringent siderite in the otherwise opaque buffer materials. As little as 0.5 percent 
siderite in the buffer could be detected and identified by the high birefringence and by the "twinkle" 
caused by the great change in relief on rotation of the microscope stage. 


SYNTHESIS OF SIDERITE 

Experimental Methods 

Synthetic siderite was used as the starting material in this study. Natural siderites were con- 
sidered undesirable because of the common partial oxidation of the material and because of the 
amounts of other divalent cations generally present in solid solution. 

Several investigators have previously synthesized siderite by combining compounds contain- 
ing Fe 2+ and C0 3 2 ' in solution, generally under high C0 2 pressure to prevent dissociation (Palache 
et al., 1944; Sharp, 1960; Graf, 1961; Powell, 1965; Johannes, 1968, 1969). Commonly used com- 
pounds are FeCl 2 , FeS0 4 , CaC0 3 , and NaHC0 3 . A second method of synthesis involving decompo- 
sition of ferrous oxalate dihydrate (FeC 2 0 4 ‘2H 2 0) has been independently developed by several 
workers (Rosenberg, 1960; French, 1964a; Weidner, 1968) after siderite was reported as an unex- 
pected product of mixes containing ferrous oxalate (Turnock, 1960). 
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Plate 1A— Synthetic siderite formed by decomposition of ferrous oxalate dihydrate in a sealed silver tube (P H 0 = 2000 bars; 374°C; 
164 hrs). Plane polarized light; largest grains are about 10/jn across. Many crystals are well-formed rhombonedra and show no trace 
of oxides. 



Plate 1 B— Recrystallization and partial decomposition of siderite to hematite and magnetite above the equilibrium temperature (Run 
125; P F = 1000 bars; 418 <, C; HM buffer; log f 02 = -21.5; 291 hrs). Plane polarized light; largest grains are 20-30 /jti across. The 
siderite has recrystallized to larger grains, but numerous grains (e.g., upper and lower right) show incipient decomposition to opaques, 
a red color, and definite magnetism. 


The siderite used in this study was synthesized by decomposition of ferrous oxalate dihydrate. 
In the synthesis runs, from 100 to 400 mg of reagent FeC 2 0 4 -2H 2 0 was placed in a 4.5-mm O.D. 
silver tube from 2 to 8 cm long. The tube was welded shut, placed in a Tuttle cold-seal bomb, and 
heated to between 350° and 380°C under 2 kb water pressure. Reaction times ranged between 6 and 
452 hours and were generally between 40 and 115 hours. Neither reaction time nor temperature 
was particularly critical to the reaction, and complete conversion of the oxalate to siderite was 
observed in virtually all runs. In shorter rims near 350°C some unconverted oxalate occasionally . 
remained, and in runs near 400°C the siderite was accompanied by small amounts of magnetite. 

The siderite used in the stability runs was entirely pure and was generally produced at about 380°C 
in runs of 90-120 hours duration. 

The siderite produced was a gray-white crystalline powder which gradually assumed a light 
buff or brownish color on drying at room temperature. The material was identified as siderite by 
both optical and X-ray methods. The individual siderite grains are colorless to pale yellow, rang- 
ing in size from less than l M m to 10 /im (Plate 1A); many of the large crystals are well- developed 
rhombohedra. This material, which showed no traces of oxides or other impurities, was used as r 
starting material for the stability rims. The value of n o , determined on six samples, was 1.875 
± 0.005, in good agreement with other determinations (Winchell and Winchell, 1951; Rosenberg, 
1960). 

Samples of the ferrous oxalate dihydrate and synthetic siderite were analyzed through the kind- 
ness of Dr. B. F. Jones and Miss S. Rettig of the U. S. Geological Survey. The partial analyses 
are (in weight percent): 



Oxalate 

Synthetic 

Siderite 

(S-31) 

Ideal 

Siderite 

Fe 

31.25 

47.5 

48.2 

Mn 

0.00 

0.05 


Na 

<0.1 

<0.1 

m 

K 

<0.1 

<0.1 


Ca 

— 

0.0 


Mg 

0.00 

0.0 


c 2 o 4 

47.2 

— 


Total 

78.45 

47.6 



These analyses show that divalent cations other than iron are not present in significant amounts in 
either the oxalate or the synthetic siderite. 
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Measurements of the d-spacmgs of the strong siderite (104) reflection indicated no perceptible 
difference in the cell dimensions between several synthetic siderites or between three siderites 
used in stability runs at varying temperatures and f 0 values Measurements were made on an 
X-ray diffractometer using a CdF 2 internal standard (Haendler and Bernard, 1951, Harker and 
Tuttle, 1955a, 1955b, Rosenberg, 1960). Measurement of the 20 1O4 value can be made to ±0.02°. 

The values of d 104 thus determined are consistent to ±0 001 A and are believed accurate to this 
value 

It was found necessary to redetermine the value of the (111) reflection of CdF 2 used as a stand- 
ard for measurements of 2 0 1O4 The original determination by film methods of the CdF 2 cell pa- 
rameters (Haendler and Bernard, 1951) gave a o = 5 3880 ± 0 0005 A., d iu = 3.103 A (20 m 
= 28 770° for Cu K a radiation). The original A.S.T.M card based on this data (5-0567) gives 
d ul = 3.1051 A (20 1U = 28 750°). Using a standard of Lake Taxaway, N. C. quartz, for which 
20 2U = 26.664° (Lipson and Wilson, 1941, Swanson et al , 1954), a value of 20 m for CdF 2 was de- 
termined as 28.701 ± 0.005° This value, which is used in the present study, is m good-agreement 
with a similar determination by Rosenberg (1960), in which a value of 20 m = 28.699 ± 0.005° was 
determined. Using additional CdF 2 reflections, a new value of a o = 5 3872 ± 0 0009 A was ob- 
tained This new value agrees, within the uncertainties, with the original value of 5.3880 ± 0 0005 
(Haendler and Bernard, 1951), but the differences in the two values of 26 nl are greater than the 
error of measurement Accordingly, the value of 20 m = 28 701° was used as a standard for meas- 
urement of the values of d 104 of the synthetic siderite 

The average value of d 104 determined from four synthetic siderites is 2.7919 A, m excellent 
agreement with previously calculated (2.7912, Graf, 1961) and measured (2.789, Sharp, 1960) values. 
A similar constancy was observed m d 104 of three synthetic siderites exposed to various tempera- 
tures and f 0 values during stability runs (see Tables 2 and 3) 2.7897 A (Run 29, HM buffer, 230°C, 

log f 0 = -34 8), 2 7913 A (Run 52, HM buffer, 318°C, log f 0 = -27.5), 2 7908 (Run 46, MGr buffer, 
493°C, log f 02 = -23.4). The latter sample was about 20 percent decomposed to magnetite. The 
average value for the three runs, 2 7906 A, is not significantly different from that of the synthetic 
siderites. 

Similar agreement was found for d-spacmgs of other reflections of the synthetic siderites 
(French, 1964a), irrespective of whether they had been used m stability runs The data indicate 
that neither variations in synthesis conditions nor changes in temperature and t 0 ^ during later ex- 
perimental studies have any effect on siderite cell size. The X-ray data further provide no evidence 
for any significant amount of Fe 3+ in siderite as a function of f 0 . 

Production of Organic Compounds During Siderite Synthesis 

A sharp, fetid odor, suggestive of organic compounds, was observed during the siderite syn- 
thesis runs. The odor often appeared upon opening the bombs, even before the sealed capsules 
containing the siderite were opened The capsules themselves were puffed, contained considerable 
positive pressure at room conditions, and produced a similar "or game" odor when opened A brief 
study of the gas phase coexisting with siderite during synthesis was undertaken to see what, if any, 
organic compounds might be forming during the process of siderite synthesis. 
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The decomposition of feri’ous oxalate dihydrate to siderite may be written 


FeC 2 0 4 2H 2 0 = FeC0 3 + CO + 2H 2 0 ( 15 ) 

In a gas phase consisting initially of H 2 0 and CO, equilibria will be established between the various 
gas species H 2 0, H 2 , C0 2 , CO, CH 4 , and 0 2 (French and Eugster, 1965, French, 1966) The suc- 
cessful production of siderite indicates that < 10 “ 2 5 bars (French, 1965, French and Eugster, 1965) 
and, if equilibrium is established with the ambient f 0 of the hydrothermal bomb (approximately 

that of the Ni-NiO buffer), f 0 ~ 10“ 30 bars (Eugster and Wones, 1962, p. 95). 

2 

Under such reduced f 0 , it is likely that CH 4 will be a significant component of the gas phase, 
together with other reduced species (French, 1966). Furthermore, an atmosphere in which H 2 0/C0 
= 2 might be expected to generate additional organic species through inorganic equilibrium reactions 
(Dayhoff et al., 1964) 

Accordingly, a brief search was made for methane and other organic species in the gas phase 
coexisting with, siderite during synthesis Dr. T. C Hoermg, of the Carnegie Institution of Wash- 
ington Geophysical Laboratory, kindly analyzed a portion of the vapor from a siderite synthesis run 
with a mass spectrometer The sealed tube was opened in the spectrometer system. Surprisingly, 
the vapor consisted almost entirely of water and complex oxygenated hydrocarbons; methane and 
C0 2 were not present in more than trace amounts. A cursory analysis of the complex mass- 
spectrometer patterns suggested the presence of propionic acid (CH 3 CH 2 COOH) and ethyl alcohol 
(CH 3 CH 2 OH), possibly present in addition were butyric acid (CH 3 CH 2 CH 2 COOH), acetic acid 
(CH 3 COOH), and ketones. These organic compounds are undoubtedly responsible for the distinctive 
odor accompanying siderite synthesis. 

The apparent synthesis of complex organic compounds in a gas phase in equilibrium with sid- 
erite has significant implications for such problems as the nature and evolution of planetary atmos- 
pheres and the prebiological formation of organic materials in meteorites and on the primitive 
earth. Numerous investigators have shown that a wide variety of organic compounds, including 
amino acids, can be produced inorganically by the action of diverse energy sources on low-pressure 
mixtures of simple gases such as CH 4 , NH 3 , H 2 0, C0 2 , CO, etc. (Miller, 1955, 1957, Miller and 
Urey, 1959, Oro, 1965, Abelson, 1966, Ponnamperuma and Gabel, 1970). Reactions of this nature 
are considered to have been necessary for producing, from chemically simple primordial atmos- 
pheres, the critical organic compounds necessary for the development of Me. 

The related problem of whether the organic compounds in meteorites have originated through 
biological activity or from abiotic synthesis has been strongly debated (see, e.g., Urey, 1966, Urey 
and Lewis, 1966; Studier et al , 1965). There is, however, considerable theoretical and experi- 
mental evidence to indicate that such compounds may form abiotically through near-equilibrium 
reactions between simple components such as H a O, CO <2 , H 2 , CH 4 , etc., in either a primordial 
solar nebula or a primitive planetary atmosphere' (Dayhoff et al., 1964, Eck et al , 1966, Studier 
et al , 1968, Hayatsu et al , 1968). Formation of organic compounds under these conditions takes 
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place by Fischer-Tropsch reactions (Oro, 1965, Studier et al , 1968) m which catalysis by iron- 
bearing phases such as Fe or Fe 3 0 4 has a definite promoting effect (Storch et al , 1951, p 1-35, 
222-308, Anderson, 1956). 

The production of significant amounts of organic compounds in a high-pressure gas phase m 
equilibrium with siderite supports the conclusion that organic compounds can be produced abiotically 
by Fischer-Tropsch reactions The results further suggest that such syntheses can occur at high 
total gas pressure corresponding to moderate depths in a planetary crust where the existence of 
CH 4 -rich phases may be favored by low f values controlled by equilibrium with graphite and 
other mmerals (French^ 1966) Siderite, m particular, may promote such synthesis catalytically, 
as the high-pressure production of hydrocarbons from reaction of H 2 with carbonate minerals 
strongly suggests (Giardmi and Salotti, 1969) 

These brief results suggest that organic compounds may form abiotically by reaction between 
a. high-pressure gas phase and siderite, other carbonates, or other iron-hearmg minerals On a 
planetary body, such suitable conditions might be attained by hydrothermal activity at moderate 
depths within the crust The organic compounds produced at depth could then be emitted to the 
surface, providing material for further biological development even in cases where such synthetic 
reactions were not possible in the planetary atmospheres or oceans themselves 


EXPERIMENTAL RESULTS 

The Equilibrium Siderite 4- Hematite + Magnetite + Gas (SHMG) 

V 

The uni variant equilibrium siderite + hematite + magnetite + gas (SHMG) is the intersection, 
of three divariant surfaces in P F - f 0 - T space, the surface, siderite + hematite + gas (SHG), 
the surface, siderite + magnetite + gas (SMG); and the degenerate buffer surface, hematite 
+ magnetite + gas (HM) The location of the SHMG curve was determined at 500 and 1000 b P F by 
reversing reactions m samples of siderite surrounded by hematite -magnetite buffers (Table 2, 
Figure 4) The values obtained are 

P F T(±10°C) - log f 0 (±0 8) 

500 363 24 7 

1000 365 24 6 


The greatest difficulty in determining equilibrium temperatures along the hematite-magnetite 
buffer was caused by the slight amount of reaction in both directions, under 5 percent. Determina- 
tion of phases was necessarily made by optical study and by the detection of color changes and 
magnetism in the sample 

« - 

Limits on. the equilibrium temperature of the isobaric invariant point SHMG at 500 bars are 
fixed by Runs 108 (354°C) and 93 (372°C) An equilibrium temperature of 363 ± 10°C is estimated 
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Table 2 


Experimental Data for Determination of the Equilibrium 
Sidente + Hematite + Magnetite + Gas (SHMG) Along the Hematite-Magnetite Buffer 


Run No 

Sample 

n 

-iogfo 

Time 

Products 

(bars) 

(hrs ) 

Sample 

Buffer 

P F = P C q 2 + P co = 2000 bars 

25 

S 

203 

37 5 

185 

S 

H + M + S 

29 

S 

230 

34 8 

328 

S 

H +M + S 

33 

S 

249 

33 0 

501 

S 

H + M + S 

30 

S 

270 

31.2 

329 

s 

H + M 

34 

S 

274 

30 8 

506 

s 

H +M+ S 

42 

S 

276 

30 7 

473 

S + H + M 

H + M + S 

38 

S 

283 

30 1 

354 

S+H + M 

H + M + (S) 

43 

S 

289 

29 6 

473 

S 

H + M 

54 

S 

290 

29.5 

710 

S + H + M 

H + M + (S) 

39 

S 

291 

29 5 

355 

S 

H + M 

53 

S 

292 

29 4 

432 

' s 

H + M 

18 

S 

301 

28 8 

91 

S + H + m 

H + M + (S) 

26 

s 

302 

28 7 

185 

S + H + M 

H + M 

131 

s 

302 

28 7 

335 

S + H + m 

H + M + (S*) 

52 

s 

318 

27 5 

415 

S 

H+ M + (S) 

55 

s 

320 

27 4 

710 

S + H + M 

H + M 

132 

HM 

254 

32 6 

335 

H + M + (S) 

H + M + (S) 

133 

HM 

306 

28 4 

331 

H + M + (S) 

H + M 


- 


P F = P cc , 2 + P co = 1000 bars 



44 

S 

241 

33 9 

238 

S 

H + M + S 

45 

s 

280 

30 5 

190 

S 

H + M + S 

58 

s 

294 

29 4 

356 

S 

H + M + S 

70 

s 

303 

28 6 

356 

S 

H + M + S 

- 71 

s 

325 

27 1 

356 

S 

H + M + (S) 

79 

s 

348 

25 6 

567 

s 

H + M + S 

49 

s 

360 

24.8 

^ 531 

s 

H + M + S 

120 

s 

360 

24 8 

372 

S + h + M 

H + M + S 

65 

s 

370 

24 2 

373 

S + H + M 

H + M + (S) 

124 

s 

380 

23 6 

* 385 

S + h + (M) 

H + M + (S) 

121 

s 

389 

23 1 

371 

S + H + M 

H + M + (S) 

59 

s 
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S + (M) 
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s 
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26.1 
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S 

H + M 
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s 
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S 

H + M + S 

93 

s 
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s 
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s 
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105 

s 
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H + M + (S) 

99 

s 
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Figure 4— Plot of Pco^ + ^co against T, showing location of the ex- 
perimentally determined univariant curve siderite + hematite +magnetite 
+ gas (SHMG). Dashed portions of the curve are extrapolated from runs 
made at 500 and 1000 bars. In each run symbol, upper box denotes be- 
havior of the siderite sample, lower box denotes behavior of the buffer. 
Displacement of some run symbols vertically is for visibility only and does 
not indicate differences in run pressure. Runs near zero pressure were 
made at 30 psl C0 2 pressure and do not locate a definite decomposition 
temperature. 
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on the_basis of these two runs The equilibrium temperature at P F = 1000 bars is bracketed by Runs 
49 (36Q°C) and 65 (370°C), the estimated equilibrium temperature is 365 ± 10°C, mdicatmg that the 
SHMG curve is virtually vertical between 500 and 1000 bars. 

The equilibrium temperature for the SHMG curve at P F = 2000 bars was not definitely deter- 
mined Definite decomposition of the HM buffer to sidente was identified at temperatures below 
280°C (Table 2), but decomposition of siderite above this temperature could not be clearly established 
Reactions were extremely sluggish-. at 2000 bars, several runs held for five weeks at about 300°C 
showed no detectable change m either direction In several runs about 375°C, the HM buffer was 
completely reduced to magnetite during the experiment, while the siderite sample showed slight, hut 
not significant, decomposition 

The reason for the anomalous behavior of the runs at 2000 bars is not definitely known The 
combination of relatively low reaction rates with high pressures may retard reaction by inhibiting 
diffusion of gaseous reactants and products through the solid sample and buffer. Such a kinetic ex- 
planation for the lack of reaction at 2000 bars is favored by the observed greater degree of reaction 
m both direction at lower pressures. Decomposition of the HM buffer to siderite appeared particu- 
larly favored by lower total pressures. 

Small amounts of siderite, 'estimated generally at about 0.1 percent, were observed rarely in 
runs made above the equilibrium temperatures at 500 and 1000 bars, m which the siderite samples 
themselves showed definite decomposition The small amounts of siderite m the buffers in these 
runs are not considered to represent stable formation of siderite at the run temperature Such 
siderite could have formed (1) during run-up of the bomb at temperatures within the stability field 
of siderite; (2) by contamination of the buffer by mechanical leakage of siderite from the sample 
In these runs, decomposition of the siderite sample was regarded as the true indicator of reaction 
direction 

Siderite may, however, also develop metastably in the buffer Since buffering results from dif- 
fusion of the fluid phase through the buffer, the portion of the buffer at the extreme ends of the tube 
will be exposed to the f 0 ^ of the bomb, which lies below that of the hematite-magnetite buffer 

(Eugster and Wones, 1962). At such lower f fl values, siderite is stable at higher temperatures 

2 

(Figure 1) and could thus form in the buffer at temperatures above the true equilibrium tempera- 
ture- of the SHMG reaction hi all the experiments, however, nothing was observed to suggest that 
such formation of siderite did occur. No color boundary was observed m the buffer, nor was there 
any apparent concentration of siderite at the ends of the sample tube. 

Such metastable formation of siderite is illustrated by Run 125 (P F = 1000 bars, T = 418°C). 

In this run, the siderite sample recrystallized into larger single crystals 5 M m to 20 in size 
(Plate IB), similar crystals of siderite are scattered through the buffer. Definite decomposition of 
the sample was indicated by a general red color and by the definite magnetism of the sample 
(Plate IB). 
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The Equilibrium Siderite + Magnetite + Graphite + Gas (SMGrG) 

-\ 

The univanant equilibrium siderite + magnetite + graphite + gas (SMGrG) is' defined by the 
intersection of three di variant surfaces* siderite + magnetite + gas (SMG), graphite + gas (the 
graphite buffer curve), and siderite + magnetite + graphite + 0 2 (SMGrG*). The latter surface lies 
entirely within the condensed region below the graphite buffer (French and Eugster, 1965) and is 
only accessible at the univanant curve Location of the univariant SMGrG curve was determined 
by the same method, samples of siderite were surrounded by a magnetite + graphite mixture, al- 
lowing reaction in both directions in a single run (Table 3; Figure 5). 

Data obtamed for the SMGr G curve are 

P F (bars) T(ilO°C) 

500 455 

1000 ' 458 

2000 * - 465 

Minor amounts of siderite decomposition were observed m all runs, even those m which sig- 
nificant amounts of siderite formed m the, buffer This decomposition, which did not exceed 5 per- 
cent, produced a dark gray color and definite magnetism in the sample. The amount of such decom- 
position appeared relatively constant with temperature and may have been produced by incomplete 
buffering during the early part of the run. 

The direction of reaction in any run was determined from (1) the absence of siderite m the 
buffer; (2) a sudden shaip increase in the amount of sample decomposition as the run temperatures 
were increased An approximate- indicator of siderite instability was the appearance of magnetite 
peaks in the X-ray pattern of the sample Above the equilibrium temperatures, the amount of 
sample decomposition mci eases rapidly (Table 3, Figure 6). Below the equilibrium temperatures 
buffer decompositions of 2-10 percent were observed With one unexplained exception (Run 85), not 
even trace amounts of siderite were observed in the buffer above the equilibrium temperature. 

At P F = 500 bars, Run 89 (448°C) shows very slight sample decomposition combined with de- 
finite development of siderite m the buffer By contrast. Run 100 (464°C) exhibits no siderite m 
the buffer, while the siderite sample was almost entirely decomposed The equilibrium tempera- 
ture lies between these runs, and the almost complete sample decomposition in Run 100 suggests 
that it lies closer to the lower temperature The equilibrium temperature estimated at P F = 500 
bars is 455 ± 10°C. 

The same criteria indicate that the equilibrium temperature^at P F = 1000 'bars is bracketed 
by Runs 130 (452°C) and 122 (466°C), the estimated value is 458° ± 10°C. 

At P F = 2000 bars, the equilibrium temperature lies between Runs 92 (462°C) and 76 (469°C) 
and is estimated as 465 ± 10°C Above this temperature, siderite samples were mcreasmgly de- 
composed, no siderite remained m a sample run at 537°C. 


- io S f o 2 ( ±0 8 ) 
25.8 
25 2 
24.4 
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Table 3 


Experimental Data for Determination of the Equilibrium 
Siderite + Magnetite + Graphite + Gas (SMGrG) Along the Graphite Buffer 


Run No 

Sample 


~ lo §> f O- 
(bars) 

Tune 
(hrs ) 

Products 

Sample 

Percent 

Decomposition 

Buffer 




P F = 

P C o 2 + p co = 2000 bars 



68 

S 

197 

41 0 

429 

S + M + g 

2-5 

M + G + S 

62 

S 

304 

32 1 

353 

S+M + g 

5 

M + G + S 

63 

s 

348 

29 6 

353 

S + M + g 

2-4 

M + G + (S) 

56 

s 

402 

27 1 

355 

S + M + g 

2-5 

M + G + S 

40 

s 

422 

26 2 

263 

S + M + g 

5 

M+G + S 

41 

s 

441 

25 4 

263 

S + M -J- g 

2 

M + G + S 

81 

s 

451 

25 0 

371 

S + M g 

2-5 

M + G + S 

92 

s 

462 

24 6 

407 

S + M + g 

2-5 

M + G + S 

76 

s 

469 

24 2 

355 

S + M + g 

10-20 

M + G 

57 

s 

470 

24 2 

356 

S+M + g 

40-60 

M + G 

51 

s 

478 

23 9 

478 

S + M + g 

90-95 

M + G 

46 

s 

493 

23 4 

338 

S+M + g 

50 

M + G 

77 

s 

537 

21 9 

355 

M + g 

100 

M + G 

47 

s 

606 

19 9 

331 

M + G 

100 

M + G 




P F = 

P C o, + p c o = 1000 bars 



69 

s 

196 

41 9 

426 

S + M + g 

3 

M + G + S 

* 61 

s 

298 

33 0 

375 

S+M + g 

5 

M + G + S 

123 

s 

421 

26 7 

372 

S + M + g 

2-5 

M + G+S 

60 

s 

446 

25 7 

356 

S+M + g 

5 

M + G + S 

130 

s 

452 

25.5 

385 

S + M + g 

2 

M + G + S 

122 

s 

466 

24 9 

372 

S + M + g 

30-40 

M+G 

87 

s 

479 

24.4 

326 

S + M + g 

70-80 

M + G + (S) 

80 

s 

515 

23 2 

475 

M + G 

100 

M + G 

128 

MG 

419 

26 8 

385 

M + G + (S) 

0 5 

M + G + (S) 

129 

MG 

491 

24 0 

385 

M + G 

0 

M + G 




IJ 

P c o 2 + P co = 500 t> ars 



67 

s 

398 

28 2 

453 

S+M + g 

5 

M + G + S 

89 

s 

448 

26 0 

406 

S+M + g 

1 

M + G + S 

100 

s 

464 

25 4 

326 

M + S + G 

90 

M + G 

94 

s 

475 

25.0 

308 

M + S + G 

80 

M + G 

101 

s 

485 

24 6 

32e 

M+ (S) +g 

98 

M + G 

85 

s 

498 

24 1 

302 

M + S + g 1 

S0-90 

M + G + (S*> 

95 

s 

525 

23 2 

307 

M + g 

100 

M + G 

90 

s 

548 

22 5 

405 

M + g 

100 

M + G * 

118 

MG 

439 

26 4 

325 

M + G + (S) 

0 5-1 

M + G+S 

119 

MG 

470 

25 2 

325 

M + G 

0 

M + G 
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Figure 5 — Plot of Pf Pco 2 + Pco against T, showing location of the experimentally 
determined univariant curve siderite + magnetite + graphite + gas (SMGrG). In each 
run symbo!,the upper box denotes behavior of the siderite sample, the lower box denotes 
behavior of the buffer. Slight vertical displacements of run symbols have been made for 
visibility and do not Indicate differences in run pressure. The dashed portion of the 
curve is extrapolated. The runs near zero pressure were made at30psi C0 2 pressure and 
do not indicate a definite equilibrium decomposition temperature. 


25 




ESTIMATED PER CENT MAGNETITE 


O SIDERITE STABLE 


) 




Figure 6— Plot of estimated percentage of siderite sample decomposition as a function of 
temperature for runs used to locate the univariant curve siderite + magnetite + graphite 
+ gas (SMGrG). The equilibrium temperature (455-465°C) are located by the sharp increase 
in decomposition with increasing temperature. Decomposition below this temperature is 
relatively constant and does not exceed 5 percent, but decomposition ?s virtually complete 
in several runs located above the equilibrium temperature and below 500°C 
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There was no observable effect of total pressure on the amount of reaction along the SMGrG 
curve, in general, virtually complete decomposition of siderite was observed m all samples run at 
30° to 50°C above tne equilibrium temperature at all pressures (Table 3; Figure 6) 

Runs made with a magnetite-graphite mixture in the sample position showed traces of siderite 
developmg below the estimated equilibrium temperatures No siderite was observed to form in 
these runs above the equilibrium temperatures, and these runs are thus consistent with the equili- 
brium temperatures deter mined from runs using siderite samples 

The temperatures obtained for the SMGrG curve in this study are supported by independent 
studies of siderite decomposition carried out by Rosenberg (1963a). In his experiments, siderite 
was decomposed under pressure m sealed gold tubes, buffering at f 0 values at or close to those of 
the graphite buffer was apparently established by precipitation of graphite or amorphous carbon 
from the gas 'phase in the tube during the run (French and Rosenberg, 1965). The temperatures at 
which a sudden increase m production of magnetite from siderite was observed are virtually iden- 
tical to those determmed by the open-tube method used m this study (Figure 7) (French and Rosen- 
berg, 1965). This agreement between results obtained by two different experimental methods 
strengthens the conclusion that the temperatures obtained do correspond to the stable equilibrium 
temperatures for the isobaric invariant points for the equilibrium siderite + magnetite + graphite 
+ gas 


DISCUSSION OF EXPERIMENTAL RESULTS 


Buffering and Metastability in a CO 2 + CO Gas Phase 

v 

It was not possible in the present study to measure directly the C0 2 /CO ratio m the gas phase 
coexistmg with the sample during an experiment, and indirect evidence must be used to indicate 
whether the buffer assemblage does m fact control the f 0 value imposed on the sample in a C0 2 
-b CO atmosphere In earlier experiments m hydrothermal (H 2 0 + H 2 ) atmospheres, buffering was 
carried out by osmotic diffusion of hydrogen through sealed platmum capsules (Eugster, 1959, Shaw, 
1963) The success of this method has been demonstrated by the reversibility of reactions and by 
the consistently varying compositions of iron-bearing solid solutions synthesized at various f 0 
values (Turnock and Eugster, 1962, Eugster and Wones, 1962, Wones and Eugster, 1965). The suc- 
cessful operation of the buffering mechanism has also been verified by direct measurement of f 
values mside the sample tube (Shaw, 1963, 1968). 

In contrast, buffering m a C0 2 + CO atmosphere such as that used in the present study involves 
mechanical diffusion of all gas species through the solid buffer and the open tubes and thus requires 
that the necessary solid-gas equilibration reactions take place in times that are short compared 
to run duration If these conditions are not met at the relatively low reaction temperatures in- 
volved (250-500°C), then equilibrium will not be attained and the determinations will not correspond 
to equilibrium siderite decomposition temperatures 
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SEALED TUBES, UNBUFFERED 
O Sid + Mt (*=5%) 

^ Sid + Mt (>20%) 


OPEN TUBES, BUFFERED 
□ Sid + Mt (-'5%), Sid FORMED 

(j Sid + Mt (>10%) 



Figure 1 — Location of the untvarianf equilibrium siderlte + magnetite + graphite + gas, as 
determined by open-tube buffered runs (French, 1964a, rectangles) and by sealed-tube ex- 
periments (Rosenberg; circles) (from French and Rosenberg, 1965). The solid line indicates 
the univariant curve established from the open-tube buffered experiments, the sealed-tube « 
experiments correspond to virtually the same curve. Stippled rectangles at left indicate the 
experimental uncertainty, ±7°C in run temperature and ± 5 percent ?n total pressure, 




Successful buffering at low temperatures by the hematite-magnetite (HM) buffer is indicated by 
several observations 

(1) Reduction of the buffer entirely to magnetite occurred during long rims at temperatures 
near 400°C, indicating that the amount of reaction between the solid buffer and the gas phase under 
these conditions is sufficient to convert the original buffer to a more reduced assemblage charac- 
teristic of the lower ambient f 0 ^. It is reasonable to assume that buffering of the sample is mam- 
tamed by the solid phases until the hematite is used up. 

(2) Distinctly different equilibrium decomposition temperatures for siderite are obtained, de- 
pending on whether the hematite-magnetite (HM) or the magnetite- graphite (MGr) buffering as- 
semblages are used If the hematite- magnetite pair did not function as a buffer, the f 02 on the 
sample would have been that of the bomb, and no more than a small difference in observed decom- 
position temperatures would have been expected - Because the f 0 values of the graphite buffer are 
close to those of the pressure vessels at these temperatures (French and Eugster, 1965), similar 
decomposition temperatures would have been observed regardless of whether the magnetite-graphite 
assemblage functioned as a buffer. The large spread in decomposition temperatures for the two 
buffer assemblages, nearly 100°C, indicates that the hematite-magnetite assemblage provided a 

distinctly different f 0 value to the samples. 

2 

(3) Both hematite and magnetite are observed as decomposition products of siderite samples 
surrounded by a hematite-magnetite mixture, although the amount of reaction is slight. Had the 
buffer not functioned, only magnetite would have been expected as a decomposition product. 

The problem of reaction and buffering by the magnetite- graphite mixture is more complex. 
Considerable experimental evidence indicates that metastably low values of the C0 2 / CO ratio may 
be preserved for significant periods, even at elevated temperatures, without precipitation of graph- 
ite (Muan, 1958, Bank et al , 1961, Bank and Verdurmen, 1963), although there is little knowledge 
about reaction rates under the high temperatures and pressures of the present experiments. 

There is little information on the rates of equilibration of C0 2 + CO atmospheres when graph- 
ite is originally present m the system The reactions are sluggish and disequilibrium may persist 
for hours or even days at low temperatures At higher temperatures, equilibration is attained or 
approached more' quickly (Bradner and Urey, 1945, see also Mellor, 1924, Remy, 1956, for sum- 
manes), but catalysis may be required to promote reaction even at temperatures of 500-600°C 

The present experimental conditions, involving relatively high temperatures (400-600°C), long 
reaction times (generally two weeks), and the presence of siderite, iron oxides, and an iron alloy 
reaction vessel as potential catalytic agents, offer optimum conditions for equilibration of the gas 
phase with the graphite buffer The following observations indicate that such buffering does m fact 
occur 

(1) Either graphite or amorphous carbon is observed to precipitate from the gas phase This 
precipitation implies that metastably low C0 2 /CO values are not maintained and that reaction be- 
tween graphite and the gas phase does establish an equilibrium f 0 value The carbon is observed 
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♦ 

outside the capsule m the bomb, but also appears in sealed capsules of siderite run unaer similar 
conditions (Rosenberg, 1963a, French and Rosenberg, 1965). 

(2) Similar graphite or carbon is observed with magnetite in samples of siderite decomposed 
above the equilibrium temperatures The material was detected by dissolving the sample m HF, a 
process which removed the finely-divided siderite and magnetite and left the carbon as a residue 
However, not enough of the carbon phase could be isolated for X-ray determination of its exact 
character. 

(3) The behavior of unbuffered siderite samples exposed to the ambient bomb f 0 was slightly 

but significantly different from siderite samples enclosed m the buffer, even when the same furnace 

and bomb were used Only a few unbuffered runs were made At 2000 bars, an unbuffered sample 

was stable at 474°C, while a buffered sample decomposed at 438°C. At P F = 500 bars, by contrast, 

an unbuffered sample decomposed at 438°C, while a buffered sample was stable at 448°C These 

results cannot be considered definitive because of uncertainties in the actual f n of the bombs and 

2 

because of the progressive change in bomb f 0 with continued use (Eugster and Wones, 1962, p 95) 
The results do imply, however, that the decomposition temperatures of siderite differ when a 
graphite buffer is employed, and the unbuffered runs provide further evidence that the graphite buf- 
fer does m fact operate under the experimental conditions Because of the similarity in f 0 values 
for the bombs and for the graphite buffer (French and Eugster, 1965), large differences m behavior 
between buffered and unbuffered runs would not be expected 

It is possible that precipitation of carbon during runs will alter the walls of the reaction vessels 
after long use, either by coatmg them with carbon or by producing carbon-bearing alloys or car- 
bides m the metal. Such processes, particularly armoring the bomb by carbon, will cause the 
ambient bomb f 0 ^ to be fixed close to that of the graphite buffer itself, thus requiring only minimal 
reaction between the atmosphere and the graphite surrounding the sample to attain the specified 
f o 2 If the carbon which precipitates from the gas phase is not pure crystalline graphite, the f 0 
produced m the bomb will be slightly different from that calculated for the graphite buffer (see 
Zen, 1963, p 934) 

These experimental observations indicate that, in a C0 2 + CO atmosphere, reaction between 
the gas phase and the solid-phase buffers is sufficient, even at relatively low temperatures, to es- 
tablish f buffering m the sample. The temperatures determined for siderite decomposition are 
thus equilibrium temperatures for the reactions indicated, a conclusion strengthened by the demon- 
stration of reversibility at each equilibrium point Reversibility was best demonstrated at the 
higher temperatures of the siderite + magnetite + graphite + gas curve, and these points are m 
good agreement with results from another experimental study which involved buffering by precipi- 
tated graphite (French and Rosenberg, 1965). Finally, the success of this buffering method and of 
the graphite buffer in particular has been well demonstrated in subsequent experimental studies, 
particularly in the determmation of the rhodochrosite stability field (Huebner, 1969) and in investi- 
gations of complex metamorphic reactions involving C-H-O gas phases (Eugster and Skippen, 1968) 
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Comparison With Other Investigators 


Johannes (1968, 1969) studied the stability gelations between Fe-Mg carbonates and aqueous 
solutions under conditions in which f 0 ^ and other volatile pressures were not specifically controlled 
Because oxalic acid was often added to his system to stabilize the carbonates, the conditions of his 
experiments are probably most comparable to the siderite synthesis runs made in the present study 
In both sets of experiments, both siderite and Fe-rich carbonate showed decomposition to magnetite 
at temperatures above about 400°C No precipitation of graphite was observed either by Johannes 
or in the present syntheses, and the decomposition temperatures observed may represent equilib- 
rium breakdown at values of f 0 intermediate between the hematite-magnetite and graphite buffers 

Seguin (1968a, 1968b) carried out studies in the system Fe-C-0-S-H a 0 involving siderite de- 
composition under a variety of experimental methods and conditions. Because of the greater com- 
plexity of his system and the presence of other volatile species, it is not possible to make direct 
comparisons between his results and those of the present study, even for his experiments in which 
oxygen buffers were used. In general, his decomposition temperatures are lower than those observed 
in the present study and are likewise virtually independent of total pressure Interpretation of his 
results is complicated by a number of experimental problems Reactions were not reversed, only 
the decomposition of siderite being noted, and many of his curves are based on relatively few runs. 
Second, his use of oxygen buffering by osmotic hydrogen diffusion is uncertain, partly because the 
diffusion rates aktemperatures of 250-500°C may not be sufficient to establish equilibrium (Eugster 
and Wones, 1962, p 94, Shaw, 1968), and partly because, with a C-H-0 gas phase m the sample, 
equilibration of f will produce different values of f 0 m the sample and the surrounding buffer 
(Eugster and Skippen, 1968, p. 503-504). Finally, some of the buffer assemblages which he used 
(e g., magnetite-iron) are clearly metastable with respect to the graphite buffer. 

Even if Segum’s (1968a, 1968b) results do represent equilibrium decomposition of siderite, the 
experiments were performed under conditions where, because of the other components in the system, 
Pco 2 + Pco < Pf 311(1 the exact values of f c0 ^ , f co , and f are not known It is therefore not possible 
to compare his temperatures with those obtained m the present study under conditions where P co 
+ P co = P F and f o 2 is specified 

The experiments of Weidner (Weidner and Tuttle, 1964, Weidner, 1968) are more comparable 

with the present investigation because they were performed m an atmosphere of CO^ + CO, using 

sealed gold tubes containing siderite samples Explicit control of f „ was not attempted, and the 

2 

observed decomposition of siderite between 500 bars and 10 kb was used to locate two reactions 
(1) siderite + hematite = magnetite + gas, (2) siderite = magnetite + graphite + gas The two 
reactions should, ideally, correspond to the SHMG and SMGrG curves, respectively, determined 
here In addition, a few open-tube buffered experiments similar to those in the present study were 
made $$ 

Decomposition temperatures determined in the sealed-tube experiments (Weidner and Tuttle, 
1964, Weidner, 1968) are significantly higher than, those observed for analogous reactions m the 
present study (French, 1965; French and Rosenberg, 1965) (Figure 8). Because of slow reaction 
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Figure 8— Location of the two univariant curves siderite +~hematite + magnetite + gas (SHMG) and 
siderite + magnetite + graphite +gas (SMGrG) as determined by open-tube buffered experiments (French, 

1964a) and by seaied-tube experiments (Weidner and Tuttle, 1964; Weidner, 1968), The greatest dis- 
crepancy is for the SMGrG curve, for which Weidner's determined equilibrium temperatures below 2000 
bars are 50-90°C higher than those determined here. The break in slope of Weidner's SMGrG curve has 
been interpreted as indicating a polymorphic transition in siderite (Weidner and Tuttle, 1964), but no cor- 
responding break is evident in the SHMG curve, 

rates at lower temperatures and pressures, Weidner did not determine the SHMG curve below about 
3 kb Extrapolating his determined curve to the lower pressures attained in the present study (Fig- 
ure 8) gives approximately similar decomposition temperatures, although the slopes of the two 
curves are distinctly different 

The disagreement in the two investigations is most severe for the reaction siderite + magnetite 
+ graphite + gas .(Figure 8). The temperatures deter mined by Weidner are from 60° to 90°C higher 
than those determined here The interpretation of the break in slope of Weidner's curve at 5 kb as 
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the result of a polymorphic transition of siderite (Weidner and Tuttle, 1964) may not be correct 
His experimental points with their uncertainties can be fitted reasonably to a single smooth curve, 
and no similar break in slope is evident m his other curve (Figure 8) The exact causes of the 
large temperature differences in the two investigations are not clear. Nor is it obvious why the 
lesults of Weidner's sealed-tube experiments differ so greatly from the apparently similar sealed- 
tube experiments of Rosenberg (French and Rosenberg, 1965) 

The sealed-tube technique itself involves a number of experimental problems that complicate 
interpretation of the results 

(1) There is no direct control of f 0 , and it must be assumed that f 0 is buffered during re- 
action as the necessary phases are produced in the sample.tube This assumption is probably more 
justified for the graphite buffer, because formation of graphite in the charge is commonly observed 
(French and Rosenberg, 1965, Weidner, 1968, p 23) However, metastably low values of f 0 can be 
produced by the dissociation of siderite, since the resulting C0 2 /CO ratio will be 2 or 3 depending 
on the decomposition product. Such ratios correspond to values of f 0 ^ that are metastably low with 
respect to the graphite buffer (French and Eugster, 1965), and if these values are maintained, sider- 
ite will appear stable at temperatures above the true equilibrium temperature (see Figure 1) 

_(2) The reactions are not reversed, and only the decomposition of siderite is used to locate 
the equilibrium temperature 

(3) Reaction rates may be reduced as a result of compaction of the sample tube and charge 
during pressurization. 

(4) Decomposition of the siderite sample during run-up below the equilibrium decomposition 
temperature is necessary to produce internal gas pressure m the sample tube. The potential 
presence of decomposition products m all runs makes it more difficult to recognize the true equi- 
librium temperature by the amount of sample decomposition Further, it must be assumed that 
the external pressure is in fact equalized by internal gas pressure throughout the entire sample, 

a condition that may not exist if some external pressure is supported by the compacted sample 
itself. 

The last two effects appear to have been important m Weidner’s experiments and may account 

for much of the difference between his results and the results reported here for the SMGrG curve 

/ 

Reaction rates in Weidner’s study appear generally lower than those observed m the open-tube 
investigations, as indicated by the negligible amount of reaction on the SHMG curve below about 
3 kb. Furthermore, the sample geometry apparently influenced his results, larger amount of re- 
action products were observed to form when larger- diameter sample tubes were used (Weidner, 
1968, p. 24-25, 38). 

A crucial point m evaluating the sealed-tube experiments is the value selected by the investiga- 
tor for the amount of sample decomposition produced during run-up Observed decomposition above 
this amount indicates that the equilibrium temperature has been exceeded, and one of the disadvan- 
tages of the sealed-tube technique is that the apparent equilibrium temperature can be a function 


33 



of the limit selected. French and Rosenberg (1965) regarded more than 5 percent decomposition 
as indicating that the equilibrium temperature had been exceeded (Figure 6). Weidner (1968, p. 26) 
generally interpreted more than 10 percent magnetite as indicating stable decomposition, although, 
m some experiments (Weidner, 1968, p 39-44), the formation of about 30 percent magnetite at 
540°C was interpreted as indicating that siderite was stable. The higher cutoff limit used by Weidner, 
combined with possible reduction of reaction rates by sample compaction, would combine to produce 
apparent decomposition temperatures in excess of the true equilibrium values. 

Weidner also performed some open-tube graphite-buffered experiments that were similar in 
sample geometry to those used here (Figure 3), Their results, however, were interpreted by him 
as supporting the higher temperature values Some of these runs were made in an atmosphere- con- 
taining 30-40 percent CO, a value much higher than permitted by the graphite buffer and correspond- 
ing to a metastably low f 0 . Equilibration of such a reduced atmosphere with the graphite buffer 
requires that the gas reaction 2 CO = C0 2 + 1/2 0 2 be promoted. The reaction rate may be slow 
enough at moderate temperatures to allow metastably low f 0 values to be preserved during the 
.experiment In such an environment, siderite will exhibit metastably high decomposition 
temperatures 

Equilibration of a pure C0 2 atmosphere with a graphite buffer, the method used here and in 
some of Weidner’s buffered experiments, is a more desirable technique The composition, in equi- 
librium with the graphite buffer at these temperatures is more than 97 percent C0 2 (French and 
Eugster, 1965), and equilibration m this case can be established with only a small amount of the 
heterogeneous reaction C0 2 + C = 2 CO. Such equilibration is much more likely during the course 
of an experiment, Weidner’s buffered experiments un a C0 2 atmosphere were made at 2000 bars. 

At 540°C, the siderite sample showed extensive decomposition (about 55 percent magnetite). At 
500°C, the siderite sample remamed stable and siderite formed in the buffer Conventional inter- 
pretation of these results would place the true equilibrium temperature between these values and 
not at the value of 556°C determined from the sealed-tube experiments. The apparent stability of - 
siderite at 500°C is an unexplained anomaly, since extensive decomposition of siderite was observed 
at even lower temperatures m the present study (Figure 6). 

/■ 

Although the causes of the differences between Weidner's results and the present determina- 
tions are not definitely established, it seems likely that much of the difference arises from dif- 
ficulties m experiment and interpretation caused by the nature of the sealed-tube method itself. 
Specific effects probably include slower reaction rates resulting from sample compaction and also 
problems m interpretation caused by the necessary production of decomposition products below the 
true equilibrium temperature Both effects would tend to produce apparent equilibrium tempera- 
tures higher than the true values. 

It is believed that the open-tube buffered methods used in the present study (French, 1965, 

French and Eugster, 1965; French and Rosenberg, 1965) are superior to the sealed-tube technique 
for determining stability relations of carbonates which involve oxidation reactions The open-tube 
technique permits direct control and estimation of f 0 Further, the open tubes remove problems 
arising from sample compaction and also allow extensive equilibration and reaction by mechanical 
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diffusion of the gas phase through the charge Most importantly, using the open-tube method, 
sample decomposition below the equilibrium temperature is negligible and the true decomposition 
temperatures can be more precisely estimated. 

Accordingly, the temperatures obtained in the present study are believed to correspond to the 
true equilibrium temperatures for the reactions studied. This conclusion is strengthened by the 
agreement between the results of both open-tube and sealed- tube experiments along the graphite 
buffer (French and Rosenberg, 1965). Weidner's results for the same equilibrium apparently 
represent maximum temperatures, and more study on the problems of interpretation and kinetic 
factors in sealed-tube experiments will be needed before the differences can be fully resolved 


SIDERITE STABILITY IN P F -f Q -T SPACE 

f o 2 

Data for the two univariant equilibria siderite + hematite + magnetite + gas (SHMG) and 
siderite + magnetite -f- graphite + gas (SMGrG), deter mined in the present study, are shown in 
Table 4. The experimentally determined temperatures are higher than those calculated from ther- 
modynamic data (see Figure 1). 

The stability field qf siderite in P F -f 0 -T space can be constructed from the experimentally 
determined positions of the two univariant curves and is conveniently presented man isobaric plot of 
log f 0 against T At P F = 500 bars (Figure 9), the siderite + gas field occupies a narrow wedge 
bounded by the three isobarically univariant curves siderite + hematite + gas (SHG), siderite 
+ magnetite + gas (SMG), and the degenerate equilibrium siderite + magnetite + graphite + gas 
(SMGrG*). Only the two isobaric invariant points have been experimentally determined, curvatures 
of the univariant curves are approximate and are based on thermodynamic calculations (French, 
1964a, Hubner, 1969) 

The graphite buffer curve (GV) divides 'the siderite + gas region into two distinct parts (French 
and Eugster, 1965) Above the curve, siderite is in equilibrium with a gas composed of C0 2 and 
CO which exerts a total pressure of 500 bars. Below the curve, siderite coexists only with oxygen 
and with an imaginary inert gas which exerts a total pressure of 500 bars. 

Along the SHG curve, siderite coexists with hematite and a gas phase of variable composition, 
the stable portion of the SHG curve terminates in the invariant point siderite + hematite + magnetite 

J 

+ gas, located at 363 °C and log f 0 = -24 7 Between this point and the isobaric invariant point 
siderite + magnetite + graphite + gas at 455°C and log f 0 = -25 8, siderite coexists with magnetite 
and a gas phase of variable composition along the stable portion of the univariant curve siderite 
+ magnetite + gas (SMG). 

The point siderite + magnetite + graphite + gas at 455°C represents the highest temperature 
at which the assemblage siderite + gas is stable at P F = 500 bars. Below the graphite buffer curve, 
the stability of siderite is determined by the degenerate equilibrium curve siderite + magnetite 
+ graphite + 0 2 (SMGrG*), which is not experimentally accessible. 
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Table 4 


Temperature and Gas Composition Data for Equilibrium Points Determined on the Umvariant Curves, With Thermodynamic Data 
for the Reaction 3 Sidente + 0 2 = Magnetite + 3 C0 2 at These Points 


P F 

T°C 
(calc ) 

T°C . 
(exp.) 

- log f o 
(bars) 

f 

I co 

(bars) 

■^co 

(bars) 

log K„ (T) 
(calc ) 

AG t ° (keal) 
(calc ) 

log K 4 (T) 
(expei ) 

AG T ° (kcal ) 
(exper ) 



Siderite + Hematite 4 Magnetite 4 Gas (SHMG) 



500 

212 

363 

24 8 

3 4 x 10" 4 

535 

426 23 

-75 7 

420 58 ± 1 0 

-59 9 

1000 

223 

365 

24 7 

8 8 x 10‘ 4 

1230 

426 42 

-75 9 

421 62 

-63 1 

2000 

244 

370* 

24 4 

1 3 x 10‘ 3 

3600 

426 78 

-76 4 

422 87 

-67 3 




Siderite 4 Magnetite 4 Giaphite + Gas (SMGrG) 



500 

234 

455 

25 7 

1.3 

570 

425 60 

-84 7 

421 12 ± 1 0 

-70 3 

1000 

277 

458 

25 2 

2 1 

1310 

425 76 

-84 9 

421 95 

-73 4 

2000 

301 

465 

24 2 

3 5 

3840 

+26 07 

-85 6 

422 95 

-77 4 


^Estimated temperature based on extrapolation from lower pressures 





300 400 500 

Temperature, °C 


Figure 9 — isobanc section at Pp = 500 bars, showing the stability field of siderite + gas as a 
function of T and log fo 2 . Thestability field is based on the experimentallydetermined isobanc 
invariant points siderite + hematite + magnetite + gas (SHMG) (363°C, log f 02 = -24.8) and 
siderite + magnetite + graphite + gas (SMGrG) (455°C, log fo 2 = -25.7). Only the part of the 
stability field above the graphite buffer curve (Gr) is experimentally accessible. 


The stability relations of siderite at total pressures of 1000 bars and 2000 bars are virtually 

identical to those at 500 bars, as the polybaric projection' (Figure 10) indicates, except for the 

* 

slightly higher temperatures of the invariant points at higher pressures The chief effect of in- 
creasing total pressure is the significant shift of the graphite + gas buffer curve toward higher 
values of f Q as the result of increasing P co ^ 
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Temperature, °C 


Figure 10— Polybaric^ro|ection of the stability field of siderlte + gas determined at 500, 
1000, and 2000 bars Pco 2 * Pco* The SHMG point at 2000 bars is extrapolated from data at 
lower pressures. The stability field of siderlte + gas is bounded by stable decompositions to 
hematite + gas, to magnetite + gas, and to magnetite + graphite + gas. The size of the field 
of siderite + gas is only slightly affected by-changing total pressure between 500 and 2000 
bars, most of the effect arises from the change in location of the graphite buffer curve (Gr) 
with changing P C02 + Pco- 


The section at P- = 2000 bars is incomplete. The SMGrG pomt lies at 465°C and log f Q = 

2 

-24.4. The SHMG point was not determined, but extrapolation of the SHMG curve from lower pres- 
sures suggests that the equilibrium temperature at 2000 bars is approximately 3 70°C, winch cor- 
responds to log f 0 = -24.2 The isobaric invariant SMG curve between these points is nearly 
horizontal 
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The stability field of siderite + gas, plotted in three dimensions in P F - log f 0 - T space 
(Figure 11), occupies a narrow wedge bounded by the three divariant reaction surfaces siderite 
+ hematite + gas (SHG), siderite + magnetite + gas (SMG), and siderite + magnetite + graphite + 0 2 
(SMGrG*). 

The graphite buffer surface (GrG) forms the lower boundary of the experimentally accessible 
portion of the field of siderite + gas. Reactions involving siderite and C0 2 become metastable at 
values of log f 0 below this surface. For this reason, the stability field of siderite + gas shown in 
Figure 11 is much smaller than that calculated m similar diagrams (Holland, 1959, Garrels and 
Christ, 1965), m which reactions involving graphite are not considered. 

A schematic isobaric section is shown in Figure 12; all the stable assemblages coexist with 
gas Above the graphite + gas curve, solid phases coexist with a gas phase composed of C0 2 and 
CO. Below the curve, the solid phases coexist with 0 2 and with an imaginary inert gas that is 
added to exert the specified total pressure The presence of the additional (imaginary) component 
allows divariant assemblages below the curve to have one more solid phase than do divariant as- 
semblages above the curve The point siderite + magnetite + graphite + gas (SMGrG) represents 
the maximum temperature of stable existence of the assemblage siderite + gas at this value of P F . 

The assemblage siderite + gas is stable over a relatively narrow range of f 0 values, below 
2000 bars, siderite is not stable above 10' 24 bars f 0 ^, and most of the siderite + gas field occupies 
an area between 10' 24 and 10' 30 bars. The quartz -fayalite - magnetite buffer curve, which is the 
upper limit of stability of fayalite (Fe 2 Si0 4 ), lies entirely in the condensed region below the graph- 
ite buffer curve over the entire temperature range for which siderite is stable (Eugster and Wones, 
1962, French and Eugster, 1965) Accordingly, the stability fields of siderite + gas and of fayalite- 
+ gas do not intersect, and the possible reaction of siderite + quartz to form fayalite is not stable. 

Changes m the value of P cc , 2 appear to have only a slight effect on the stability of siderite be- 
tween 500 and 2000 bars, the equilibrium temperatures of the two umvanant curves vary less than 
10 °C in that interval. By contrast, siderite stability is strongly affected by changes in the value of 
log f Q For example, at 500 bars, a change m log f 0 of 1 0 will change the equilibrium tempera- 
ture of the assemblage siderite + magnetite + gas by more than 50°C. This effect is more striking 
at higher total pressures, where the SMG curve becomes flatter Small changes in f will produce 
similarly large temperature variations m the assemblage, siderite + hematite + gas. 

The isobaric sections demonstrate that, regardless of' total pressure, an increase mf 0 at 
constant temperature favors the decomposition of siderite to hematite or magnetite. The attitude 
of the SHG surface indicates that a temperature increase at constant P F and f 0 promotes the for- 
mation of siderite from hematite However, because of the negative slope of the SMG curve at con- 
stant P F , an increase in temperature under the same conditions will promote decomposition of sid- 
erite to magnetite. The SMG curve appears to be an exception to the general observation that an 
increase in temperature at constant f _ favors formation of the more reduced of two equivalent 
assemblages (Eugster, 1959, Eugster and Wones, 1962) Some of the analogous decomposition 
curves for rhodochrosite at constant P F have similar orientations (Huebner, 1969, Figures 8-9) 
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Figure 11 — Block diagram In Pp - log f 02 - Tspace showing the stability field of siderlte 
+ gas and its stability relations. The stability field of slderite t-gas (ruled areas on sec- 
tions) is a narrow prism bounded by di variant reaction surfaces (shown only as univariant 
curves on sections) involving decompositions to hematite + gas (SHG) r magnetite + gas 
(SMG), and magnetite + graphite + gas (SMGrG + ). The equilibrium surface siderite 
+ magnetite + gas (stippled) is bounded by the two univarianf curves determined experi- 
mentally. Stable mineral assemblages outside of the siderite stability field are indicated 
by larger lettering, smaller letters (HM, Gr) indicate buffer curves. 
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Figure 12“Schematic isobaric section through the system Fe-C-O, showing stable assemblages 
as a function of log fa 2 and T. All solid phases coexist with a gas phase. Stable univariant equi- 
libria are shown by solid lines, metastable extensions are dashed. Boxes indicate the assemblages 
stable in each region. Above the graphite buffer curve (Gr), Pp — Pco 2 + Pco/ the gas phase is 
variable m composition, and an assemblage of one solid + gas is divariant. Below the graphite 
buffer curve, Pf = Po 2 , and two solids + gas form a divariant assemblage. 


COMPARISON OF EXPERIMENTAL RESULTS WITH THERMODYNAMIC CALCULATIONS 


The experimentally determined univariant curves generally lie about 150°C above curves cal- 
culated from thermodynamic data {Table 4). The experimental curves have considerably steeper 
slopes (higher values of dP F /dT) than do their calculated counterparts This difference may reflect 
in part the uncertainty m the equilibrium temperature. When the curves are steep, small 
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uncertainties in the temperature determinations produce disproportionate uncertainties in the cal- 
culated slopes (see Burnham and Jahns, 1962, Figure 7). 

Because all the experimentally determined isobaric invariant points lie on the divariant reac- 
tion surface, sidente + magnetite + gas, the experimental data may be used to calculate values of 
AH° and AG° for the decomposition of sidente to magnetite (Equation 4) by computing values for 

logK 4 = 3 log f C o 2 “ J log fo 2 

and applying the relationships. 

AH° = - 2 303 R(d logK/d(l/T>) (16) 

AG° = - 2 303 RT log K (17) 

hi theory, such experimentally determined values of AH° and AG° would allow a more accurate de- 
termination of the values of AH£ T and AG f ° T ' for siderite Eugster and Wones (1962) applied the 
same method to calculate thermodynamic data for the iron-biotite, anmte. 

Calculated and experimentally determined values of log K 4 (Table 4) are presented in Fig- 
ure 13, together with values of the quantity (3 log f co - 1/2 log f 0 J specified by the hematite- 
magnetite and graphite buffers as a function of temperature The experimental values of log K 4 are 
quite closely grouped and are considerably smaller than values calculated from AG° for 1 bar and 
2000 bars. The lines which connect the experimental points at the same total pressure have high 
negative slopes which become more negative at lower pressures. This behavior contrasts with the 
positive slopes of the calculated lines. The negative slopes observed for the lines through the' ex- 
perimental points imply that the value of AH° is positive. Positive values of AH 0 for the experimental 
data are also suggested by the negative slope of the SMG curve in the isobaric sections (Figures 9 
and 10), since the value of (e log f 0 /d T) p 'will tend to have a sign opposite to that of AH° (French, 
1964a, p 48). 

Values of AH° derived from the experimental results are positive and range from +7 to +12 
kcal/mole between 500 and 2000 bars, in sharp contrast to the values of -13 to -16 kcal/mole de- 
rived for these total pressures from thermodynamic data (Figure 13) These differences are too 
great to be caused by the effect of changing total pressure on AH° (see Orville and Greenwood, 1965), 
because the differences exist between calculated and experimental data for the same total pressure, 
and because the effect of total pressure on this reaction can be calculated to be only about -2 keal/ 
mole between 1 bar and 2000 bars. 

Similar differences exist m values of AG° calculated from the experimental results The values 
of AG° calculated from thermodynamic data (Table 1) are about -75 kcal/mole on the hematite- 
magnetite buffer and about -85 kcal/mole on the magnetite-graphite buffer (Table 4) The experi- 
mental results are about 10 kcal/mole more positive about -64 kcal/mole (SHMG) and -74 kcal/ 
mole (SMGrG) Assuming that AS 0 is constant with temperature, the experimental results yield a 
value of AS° - 111 cal/mole-deg, in contrast to a calculated value of about +96 cal/mole-deg 
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Figure 13 — Comparison of calculated and experimentally determined values of AH 0 for the reaction* 3 si derite + 1/2 
O 2 = magnetite + 3 CO 2 {Equation 4), derived from plotting the equilibrium constant against 1/T. Values of log 
(T) ( s log Kb) derived from thermodynamic data are shown by solid lines for P F = i bar and 2000 bars. Dashed 
lines indicate values of the quantity 3 log f CO 2 - 1/2 log fo 2 as a function of temperature, using fo 2 values along 
the hematite-magnetite (HM) (long dashes) and graphite (Gr) (short dashes) buffers. Intersections of the dashed lines 
with the solid lines correspond to calculated equilibrium temperatures for a particular pressure. 

Values of log K 4 (= log K e ) calculated from experimental results (Table 4) are plotted as solid circles (the open 
circle indicates an extrapolated value for the SHMG point at 2000 bars). Boxes around the points indicate the esti- 
mated uncertainty (log K 4 ±1.0; 1/T x 10 3 = 0.03). Values of AH 0 obtained from the experimental equilibrium con- 
stants by connecting points at the same, total pressures differ significantly from the calculated values and are opposite 
in sign. However, the fact that the error boxes overlap the entire range of variation in log K 4 suggests that the de- 
rived values of AH 0 may not be meaningful. * 
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Using the relationship 


( 18 ) 


AG 2 ° 98 o k = AG t ° + AS 0 (T-298) 

O 

the experimental data yield a value of AG 29gOK = -25 5 kcal/mole and AH 0 = +7 6 keal/ mole, m con- 
trast to the calculated values of -42 9 kcal/mole and -13 9 kcal/mole, respectively (Table 5) 

Weidner (1968) did not calculate thermodynamic data from his experimental results. Approxi- 
mate values of AG° and AH° can be calculated from his data by assuming that f co ^ = P F (French 
and Eugster, 1965) and neglecting the change in f 0 of the buffers with changing P F (Eugster and Wones, 
1962) The errors introduced by these two effects are opposite in sign and will probably not pro- 
duce uncertainties greater Than those m the experimental data themselves. With these assumptions, 
Weidner's (1968) data yield values of aG° 298 = -29 1 kcal/mole, AH° = -1.4 kcal/mole, AS° = +93 

cal/mole-deg (Table 6). 

* 

Values of the thermodynamic quantities for the decomposition of siderite to magnetite (Equa- 
tion 4) derived from these two experimental studies are given in Table 6, together with values cal- 
culated by other workers (Holland, 1965, Yui, 1966) for the same decomposition The latter two 
values are not completely independent, having been based on the earlier results of French and 
Rosenberg (1965). The results calculated from experimental data appear to indicate a higher value 
of ASjgg and correspondingly more positive values of AG° 298 and AH ° 298 than are obtained directly 
from the thermodynamic data themselves. 

Because of the large effect of experimental uncertainties on the calculation of the thermody- 
namic quantities from the experimental results, it is doubtful whether the differences obtained are 
meaningful or whether these calculations can be used to evaluate the experimental data themselves 
The operation of calculating AH° (Equation 16) involves the quotient of two small difference terms 
and even reasonable estimates of experimental uncertainty (T eq ± 10°C, log K 4 ± 1.0) produce un- 
certainties as large as several hundred percent in the derived value of AH°. 

The corresponding values of AG° are less affected and the stated experimental uncertainties 
produce an uncertainty m AG° of about 5-8 percent or about ±5 kcal/mole, comparable to the un- 
certainty m the value of AG° derived from ther- 
modynamic data (Table 1). Even this moderate 
undertamty can produce large differences in 
calculated equilibrium temperatures (see Fig- 
ure 2) However, the subtractions and other 
operations necessary to derive values for AG ° 

, 2 9 8 

(Equation 18) multiply this uncertainty to the 
point where the calculated value of AG° 9S has 
a minimum uncertainty of about 100 percent. 

Because of the large uncertainties, the dif- 
ferences between the values in Table 5 is prob- 
ably not significant. The difficulty of using the 
experimental data to evaluate the experiments 
themselves is indicated further by the fact that 


Table 5 


Calculated and Experimentally Determined 
Values of AH 0 for the Reaction 
Siderite + 0 2 = Magnetite + 3 C0 2 


P F (bars) 

AH° 

(keal ) 

Calculated 

Experimental 

1 

' -13 9 ± 7 0 

— 

500 

— 

+12 5 

1000 

— 

+7 6 

2000 

-16 1 ± 7 0 

+1 8* 


•Based on estimated temperature of 370 °C for the SHMG curve at 
2000 bars 


44 



Table 6 


Calculated Thermodynamic Data for Sidente Formation and Decomposition to Magnetite at 298°K. 


(1) Reaction 3 FeC0 3 + l/2 0 2 

= Fe 3 0 4 + 3 CO 2 



(kcal ) 

AHjgg 

(kcal ) 

as 298 

(eal/mole-deg) 

calculated (French, 1964a and 
this study) (Table 1) 

-42.87 

-13 90 

+93 0 

experimental (French, 1964a and 
this study) 

-25 49 

-7 60 

+111 1 

Weidner (1968) 

-29 10 

-1 39 

+93 0 

Holland (1965)* 

-22 12 

+7 43 

t 

+99 3 

Yui (1966)* 

-29 7 

-0 14 

+99 3 

(2) Reaction Fe + C + 3/2 O a = Fe C0 3 


AG° f 298 

(kcal/mole) 

AH f 29 8 

(kcal/mole; 

AS 2 ° 98 (siderite) 
(cal/mole-deg) 

calculated (Robie, 1966) 

-161 06 

-178 20 

23 9 

experimental (French, 1964a and 
this study) 

-166 6 

-185.4 

17 7 


‘Recalculated from author’s original values for the analogous reaction 3 FeCO - — Fe^O^ + 2 CO , + CO 


the strongly divergent experimental results of French (1964a) and Weidner (1968) give comparable 

s. 

values for the thermodynamic parameters at 298 °K 

The large uncertainties in thermodynamic values obtained from the experimental data indicate 
•that any new thermodynamic data for sidente calculated from the experimental results will probably 
not be meaningful If the data in Table 6 are, however, taken at face value, then the differences be- 
tween experimentally derived thermodynamic quantities and those calculated from original thermo- 
dynamic data arise from differences in the thermodynamic quantities of one or more of the partici- 
pating phases The uncertainty m AG f ° 298 for magnetite and CO, is probably less than ±2 kcal 
(Coughlin, 1954) If the entire difference is assigned to the ther mo dynami c. data for sidente, a 
change of about +7 kcal/mole in AH f ° 29S and of about +5 kcal/mole in AG f ° 298 at 600-700°K would 
bring the experimental and calculated values into agreement (Table 6). This would correspond to 
a change of about +4 percent in the present values (Kelley and Anderson, 1935, Robie, 1962). 

Thermodynamic data for siderite at elevated temperatures (600-700°K) are not definite enough 
to exclude a difference of this magnitude. Current heat-capacity data (Kelley, 1960) are based on 
studies between 50° and 300°K of a natural sidente containing only 90 percent FeC0 3 (Anderson, 
1934), values of C p were corrected by Anderson for the additional components CaC0 3 , MgC0 3 , and 
MnC0 3 m the sample Extrapolation of these data to higher temperatures could result m the in- 
dicated difference. 
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High-temperature heat- capacity data on pure siderite will probably be needed to resolve the 
question entirely It should be noted that estimated values of C p for rhodoehrosite (MnCO s ) based 
on Anderson's data (Kelley and Anderson, 1935) are 5 to 10 percent lower m the range 500-700°K 
than are values determined by direct measurement (Moore, 1943), and that similar differences are 
observed in values of AG f ° T for rhodoehrosite determined from experimental stability data (Huebner, 
1969). It is very likely that a similar difference could exist m the thermodynamic data for siderite. 


GEOLOGICAL APPLICATIONS 


Data on the stability of siderite may be applied to estimate conditions present during the meta- 
morphism of iron formations (French, 1968) and the emplacement of siderite -bearing hydrothermal 
veins. The experimental results determined here are in a simpler system than present in nature, 
and application of the experimental data to natural processes is limited by the additional components 
present m the natural environment. 


(1) Water is undoubtedly present m significant amounts many gas phase coexisting with natural 
siderite. Introduction of water into the gas makes possible the formation of H 2 and CH 4 (French, 
1966) and produces a situation in which P C02 + P co < P F . These effects must be considered in 
estimating the decomposition temperatures for natural siderite. 


The effect of conditions where P F / P C02 may be approximately calculated for the SHMG and 
SMGrG univariant curves from the relation 


AV, +(v ro ) 

AS 1 jJpE co 2 

V dT L ~ AV ~ dT * AS 


(19) 


(Thompson, 1955; Greenwood, 1961, p. 3924-3925). Applying this relation to the univariant curves 
gives 



( 20 ) 


Values of dP t /dT for both umvariant curves correspond to slopes of about +5 deg/1000 bars. 
Calculated values of the ratio AV/AV S are about -2 to -4. Thus the values of (<?P s /<9T) Pe are 
approximately -40 to -80 bars/deg 


>co 2 


The corresponding umvariant curves for fixed values of P Eco will lie at temperatures below 
the umvariant curves which are specified by the condition P t = P c0 ^ (Thompson, 1955, Greenwood, 
1961, Figure 1). The umvariant curves will have steep negative slopes which correspond to a change 
of -10° to -20°C for each 1000 bars increase in P s 

It appears that the decomposition temperature of siderite will be affected only slightly by 
changes in total pressure at constant P EcC) . The decomposition of siderite to magnetite thus 
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represents a geothermometer which is strongly affected by changes in f 0 but which is relatively 
insensitive to changes m both P Eco and P t . 

A second effect of water m the gas phase is the potential stabilisation of iron hydroxides such 
as goethite and lepidocrocite (FeO-OH), particularly at lower temperatures. Available stability 
data on the hydroxides (Deer et al., 1962, v. 5, p 118-127) are not sufficient to evaluate this effect 
m detail 

(2) Natural siderites commonly contain 10 to 20 mole percent'of other components, chiefly 

MgCO a and MnC0 3 , m solid solution. The decomposition curve for pure magnesite (MgC0 3 ) lies 

at 750° to 850°C for values of P co between. 500 and.2000 bars (Barker and Tuttle, 1955a) and is 

virtually independent of f Q , while the stability field of rhodochrosite (MnC0 3 ) extends to higher 

temperatures and f 0 values than does the field of siderite (Huebner, 1969). The stability relations 

of the various carbonate minerals whose components can substitute in siderite (Figure 14) indicate 

that the stability of siderite will be extended to higher temperatures and f n values by the addition - 

2 

of Mg 2+ , Mn 2+ , or Ca 2+ . The decomposition curves for various intermediate carbonate composi- 
tions will be oriented at diverse angles within the envelope established by the stability relations of 
the pure end-members (Figure 14, also see Huebner, 1969, p 477, Wanes and Eugster, 1965, p. 
1254). 

Despite these limitations, the present experimental data may be applied to natural occurrences 
of siderite. 

The stability of siderite is strongly dependent on the value of f 0 . Siderite is m stable equilib- 
rium with gas over a narrow range of f 0 and is not stable above 10" 2 4 bars f 0 below P F = P cc , 2 
+ P c0 = 2000 bars Siderite coexists stably with magnetite and gas over a very narrow interval of 
f 0 values which correspond to gas compositions in which the values of P cc , 2 /P c0 exceed 500. The 
formation of oxides from siderite at constant P co ' + P co may occur either through increasing tem- 
perature or through increasing f 0 values in the environment 

Siderite-hematite assemblages (presumably primary) are occasionally found in sedimentary 
iron formations (Gruner, 1946, p. 31, French, 1968, p. 29-30). Metamorphism of iron formations, 
however, apparently occurs at T and f values within the magnetite stability field, and thus the 
assemblage siderite + magnetite + gas assumes major importance m evaluating such rocks For 
experimental conditions, which cover the values of P co in most geological environments at moder- 
ate depth, the decomposition of siderite to magnetite may oecur between 363° and 465°C, depending 
on the values of both f 0 and P Eco m the coexisting gas phase. For pure siderite, 465°C is the 
maximum temperature of stable existence between 500 and 2000 bars P C02 + P co . 

Geological studies of metamorphosed iron formations indicate that iron-rich siderites generally 
do not decompose directly to magnetite Instead, they react with the available quartz and water to 
_ produce amphiboles rich m the grunente (Fe 7 Si 8 0 22 (OH),) end- member. Such reaction has been 
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Figure 14 Isobaric sechon at Pp = Pco 2 Pco = 2000 bars, showing stability relations of siderite and other carbonates as a 
function of log f 02 and 1/T. Heavy lines indicate carbonate decomposition reactions, light fines indicate oxygen buffer 
curves. The condensed region below the graphite buffer curve is indicated by ruling. The stability field of rhodochrosite 
(MnCC> 3 } (Rh) + gas (Huebner, 1969) is circumscribed by decompositions to several Mn oxides and extends to higher tem- 
peratures and f o 2 values than does the stability field of siderite (Sid) + gas. The vertical lines indicate decomposition tem- 
peratures for carbonates whose decomposition forms only the 1 1 oxide and are therefore independent of f 02 smithsomte 
(ZnC0 3 ) (Sm) (Harker and Hutta, 1956), magnesite (MgC0 3 ) (Ms) (Harker and Tuttle, 1955a), and calcite (CaC0 3 ) (Ct) 
(Harker and Tuttle, 1955a). The relative positions of the stability fields indicate that the introduction of Mn 2 + , Mg 2+ , and 
Ca 2+ into siderite will stabilize the resulting Fe~rich carbonate to higher temperatures and f 0 values. The effect of in- 
cluded Zn 2+ will be to stabilize the carbonate at higher f 02 values with only a minor change in stability temperatures 




observed in many regions the Lake Superior district (Van Hise and Bayley, 1897, p. 368, Irving 
and Van Base, 1892, Allen and Barrett, 1915, James, 1955, French, 1968), South Dakota (Gustafson, 
1933), India (Rao, 1934), and Britain (Tilley, 1938). 

Formation of grunerite from sidente and quartz may be written 

7FeC0 5 * 8 Sx 0 2 + H 2 0 = Fe 7 Si g 0 22 (0H) 2 + 7 C0 2 (21) 

The equilibrium depends on both i CQ ^ and f H Q , but is independent of i 0 ^ for a purely ferrous 
grunerite. On a plot of f 0 vs T, for specified i co ^ and f H 0 , the reaction would appear as a line 
parallel to the f 0 ^ axis, intersecting the siderite stability field at some temperature below that of 
maximum stability. 

Magnetite is a commonly associated mineral in such rocks, suggesting that two other related 
reactions are possible under the same conditions 

3FeC0 3 + k 0 2 = Fe 3 0 4 + 3C0 2 (4) 

7 

7 Fe 3 0 4 + 24 Si0 2 t3H 2 0 + 3 Fe, Si 8 0 22 (0H) 2 i-jOj (22) 

Formation of grunerite m magnetite-quartz assemblages (Equation 22) is observed in iron forma- 
tions at the same level of metamorphism (Miles, 1943, 1946, James, 1955, Marmo, 1956, Figure 2), 

The present experimental data on siderite stability indicate that siderite is not stable above 
465°C, on this basis, an approximate temperature of about 300° to 400°C is estimated for the forma- 
tion of grunerite from iron carbonates These temperatures correspond approximately to the equi- 
librium between siderite and magnetite under the experimental conditions, and they will vary 
slightly in the natural environment, depending on f co ^, f H D , an d the amounts of other cations present 
in the carbonate These estimates are consistent with geological observations. James (1955) ob- 
served that grunerite develops m iron formations approximately coincident with the appearance of 
garnet in associated pelitic rocks, he estimates a temperature of about 300°C for "the garnet isograd 
in this region, 

Grunerite is invariably produced by metamorphism of siderite-quartz assemblages, indicating 
that water is generally available. No occurrence of fayalite produced anhydrously from the reac- 
tion of siderite + quartz has been recognized. As mentioned earlier, the stability fields of pure 
siderite and fayalite do not intersect under the experimental conditions, and, in the absence of water, 
a siderite + quartz assemblage would be converted to magnetite + quartz, which could then react 

to form fayalite at sufficiently low values of f 0 (Yoder, 1957, Gundersen and Schwartz, 1962, French, 

. 2 

1968). These conclusions apply rigorously only to pure siderite and fayalite. Intermediate Mg-Fe 
carbonates (Huebner, 1969) and olivines (Fisher, 1967) will be stable at higher values of f 0 Reac- 
tions between intermediate Mg-Fe carbonates and quartz to form Mg-Fe olivmes may be stable at 
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f 0 values above the quartz-fayalite-magnetite buffer curve. Such reactions may explain the oc- 
casional occurrences of siderite- olivine assemblages in metamorphosed iron-rich rocks {Tilley, 
1936, Klein, 1966) 

In localities where quartz is absent from siderite-beaxing rocks, magnetite is produced by 
metamorphism. Goodwin (1962) has described magnetite-rich aureoles produced in a pure siderite 
bed by a diabase dike at Michipicoten, Ontario A minor amount of grunente is also found m the 
aureole, associated with magnetite and quartz m the non-carbonate iron formation. Temperatures 
above 400°C at this locality are likely on both geological and mmeralogical grounds. 

Hydrothermal veins commonly contain siderite as a late-stage mineral (see e.g., Londgren, 

1933, Shaw, 1959). The well-known siderite-rich vein at Roxbury, Connecticut (Silliman, 1820) con- 
tains siderite and quartz without grunerite, suggesting a temperature of formation below 300-400' > C. 
The rocks which enclose the vein have been metamorphosed to the. sillimamte grade (Gates, 1959), 
suggesting that the vein was emplaced after the main period of metamorphism. Similar temperatures 
for formation of siderite- bearing veins have been estimated from other experimental data on the 
stability of Fe-Mg carbonates themselves (Johannes, 1968, 1969) on the basis of the instability of 
siderite and Fe-rich carbonates above about 400°C in aqueous solutions 


SUMMARY AND CONCLUSIONS 

■ ,>■ 

The stability of siderite, FeCO a , depends on the value of f 0 as well as on P cc , 2 and T. ‘Deter- 
minations of siderite stability were carried out m an atmosphere of CO a + CO, using solid-phase 
buffer assemblages to control f 0 . 

The stability field of siderite + gas is circumscribed by stable decompositions to (1) hematite, ^ 
(2) magnetite, (3) magnetite + graphite The latter reaction lies m a condensed region below the 
graphite buffer curve which is not experimentally accessible Reactions of siderite to form iron „ 
and wustite are metastable. 


The univariant equilibrium curve siderite + hematite + magnetite + gas (SHMG) passes through y 
the points 363 °C at P F = 500 bars, and 365°C at P F = 1000 bars. The umvanant equilibrium ! 

curve siderite + magnetite t graphite +„ gas (SMGrG) passes through the points 455°C (500 bars), I 
458°C (1000 bars), and 465°C (2000 bars). The experimental curves lie at considerably higher tem-7 
peratures than do those calculated from available thermodynamic data.* Temperatures for the same 
equilibria derived, from sealed-tube experiments (Weidner and Tuttle, 1964, Weidner, 1968) lie at 
even higher temperatures Thermodynamic analysis of the experimental data indicates that these 
differences could result entirely from experimental uncertainties The differences could also be 
produced if the value of AG f ° T for siderite at elevated temperatures is about 4 percent more posi- 
tive than the presently accepted value. The quality of the available thermodynamic data for sid- 
erite indicates that such a difference is quite possible 

Under the experimental conditions of the present study, the stability field of sidente + gas is 
entirely separated from the stability field of fayalite, and reactions between siderite and quaitz to 
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produce fayalite are metastable. This observation agrees with the apparent absence of natural 
siderite-fayalite assemblages The addition of Si0 2 to the system Fe-C-0 thus has no effect on 
the stability field of siderite + gas beyond the addition of quartz to all assemblages 

Geological application of the experimental results is complicated by the fact that (1) natural 
siderites generally contain other elements in solid solution, and (2) natural gas phases during meta- 
morphism contain considerable H 2 0, so that P co + P co < P F Maximum stability temperatures 
for pure siderite between 500 and 2000 bars P co + P co are established by the univariant curve 
siderite + magnetite + graphite + gas, which passes through the mterval 455-465°C. The effects of 
conditions where ^ C o 2 + p co < p f "Will be to lower these temperatures slightly The- stability of 
siderite is strongly dependent upon f 0 , slight variations in this quantity will stabilize the assemblage 
siderite + magnetite + gas over a temperature interval of about 100°C 

Studies of metamorphosed iron formation indicate that siderite in these rocks reacts with 
available quartz and water to form grunerite-rich amphiboles The reaction probably occurs at 
temperatures between 300° and 400°C Coexistence of siderite and quartz (as in many hydrothermal 
veins) indicates probable temperatures of formation below these values Formation of magnetite 
dunng'contact metamorphism of quartz-free sidente-rich iron formations probably occurs at tem- 
peratures above 400°C 
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